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ABSTRACT 

Many strong gravitational lenses lie in complex environments, such as poor groups of galaxies, 
that significantly bias conclusions from lens analyses. We are undertaking a photometric survey of 
all known galaxy-mass strong lenses to characterize their environments and include them in careful 
lens modeling, and to build a large, uniform sample of galaxy groups at intermediate redshifts for 
evolutionary studies. In this paper we present wide-field photometry of the environments of twelve 
lens systems with 0.24 < zjens < 0.5. Using a red-sequence identifying technique, we find that eight of 
the twelve lenses lie in groups, and that ten group-like structures are projected along the line of sight 
towards seven of these lenses. Follow-up spectroscopy of a subset of these fields confirms these results. 
For lenses in groups, the group centroid position is consistent with the direction of the external tidal 
shear required by lens models. Lens galaxies are not all super-L^^ ellipticals; the median lens luminosity 
is ^ i*, and the distribution of lens luminosities extends 3 magnitudes below (in agreement with 
theoretical models). Only two of the lenses in groups are the brightest group galaxy, in qualitative 
agreement with theoretical predictions. As in the local Universe, the highest velocity-dispersion (cr) 
groups contain a brightest member spatially coincident with the group centroid, whereas lower-tr 
groups tend to have an offset brightest group galaxy. This suggests that higher-tr groups are more 
dynamically relaxed than lower-cr groups and that at least some evolved groups exist by z ~ 0.5. 
Subject headings: gravitational lensing — galaxies: clusters: general — quasars: general — galaxies: 
evolution — galaxies: photometry 



1. INTRODUCTION 

Strong gravitational lensing of active galactic nuclei by 
intervening galaxies is a versatile tool for studying many 
astrophysical issues. Measurements of the time delay in 
brightness fluctuations between different images, when 
combined with an assumed mass model, give a measure- 
ment of the Hubble con stant indepen dent of the canoni- 
cal distance ladder fe. g..lRefs dal 196 l:lVanderriest et al.l 

" iBlandford & Naravan' '19921 
ICourbin. Saha. & Schcchtei 



1989*; 'Roberts et a^ 
Kgeton & Ko chan gk 
2002t iKochanek fc 




Schechterl l2004t ISchechterl [20(1^ . 
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Gravitational lensing statistics can be used to constrain 
the density and the equation of state of dark energy 



e.g., rTurner|ll9 90: Kochanck 1996a; Co orav fc Huteren 
iChad 12003; .Mitchell et al. .2005.) . Strong lens- 



ing also allows the determination of the m ass, shape 
and evolution of dark-matter halos (e.g., 'Kochaney 
199 It iK ceton et al. 1998; Trcu fc Koopmans 2002a ij 
Rusin et alJ l2003t iRusin fc Kochanekl l2005|) . Anoma- 
lies in the flux ratios of individual lenses are the ba- 
sis fo r studies of subs t ructu r e in dark matter halos 
(e.g., [M ao fc Schneider 199^; 'Metcalf fc Madaul l200ll 
iDalal fc Kochanck 2002; Metcalf et al. 2004), measure- 
ments of the size of AGN accretion disk sizes (e.g.. 
Ranch fc BlandfordlTMt iWvithe et"alll2000t lYoneharal 
and the size of AGN broad-line 



regions (e.g., lAbaias et a,lJ HHHl ILewis fc Tba,tal HM 
■Richards et al.ll2004|) . 

These studies are all sensitive to the mass 
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model used to calculate the lensing potential. 
Understanding whether the lens potential is sig- 
nificantly affected by a complex environment is 
crucial for obtaining ac curate models and avoid- 
ing significant biases l)Keeton fc Zabludofj 120041) . 
Simple arguments suggest that at least 25% of 
lens galaxies should lie in bound groups or clus- 
ters of galaxies (jKeeton. Christlein. fc ZabludofH 
I2000|l . but most lens environments have not 
yet been studied. A few gravitational lenses 
are already known to lie in galaxy groups (e.g., 
Kuudic et al.l ITflflTi ITflM iTourv fc Kocha^ 

19991 iFassnacht et al.l '2006') and clusters (e.g.. 
Young et al.1 I198H Fisch er. Schade. fc BarrientosI 
19981 iKneib. Cohen. fc~Hiorthl 12000(1 . Extended X-rav 
emission, indicative of more massive, dynamically 
evolved groups has been detected around some lenses 
(PG1 115-H080 and CLASS B1422-H231; IGrant et all 
l200l . but not arou nd others (CLASS B 16 00-^434 and 
CLASS B1608-^656: lDai fc Kochalleai200l . 

Attempts have been made to correlate lens 
models wit h lens environments, w ith mixed re- 
sults (e.g.. iHogg fc BlandfordI 119941 iSchechter et~an 
19971 iKeeton fc Kochanekl 119971 iKundic et all Il997al 
Kundic et alJll997bt iLehar et alJl2000|) . However, these 
studies typically examine a region only ~ 30" (~ 130 
kpc at z = 0.3 ) around the lens, and group virial radii 
(- 500/1-1 kpc: lZabludoff fc Mulchaevlll99^ can extend 
out to a few arcminutes. Most recently, Momcheva et 
al. (2006; hereafter IMQ0) find that six of eight lenses in 
that spectroscopic survey lie in groups. Some of those 
lenses were targeted because a group was pre viousl y 
known or suspected, so it is not yet clear whether IM06I 's 
high group fraction is representative of most lenses. A 
larger, less-biased survey of lens environments is clearly 
required. 

As gravitational lensing involves an integration of mass 
along the line-of-sight, one must also consider the pos- 
sible effect of line-of-sight galaxy groups, both fore- 
ground and background. If the impact parameter is 
small and the mass large, an interloping group can 
have a significant impact (> 0.05 in the normalized 
shear and convergence terms) in lensing models (M06il . 
although theoretical estimates of whethe r this effect 
is common yield conflicting results (e.g., 'Seliak' '1994 
iBar-Kana 1996; Kccton ct al. 1997; Prcmadi fc MarteJ 
l2004» . Observational studies find evidence of line- 
of-sight groups towa rd a few gravita tional lenses, 
such as B0 712-H472 llFassnacht fc Lubi n 2002), MG 
J113H-0445 fTonrv fc Kochane kll200(f ~and B1608-H656 
(p^assnacht et al. 2006). In M06 we present spectroscopic 
identification of a massive group behind MG J0751+2716 
that significantly impacts that lens model. A detailed 
census of line-of-sight groups is clearly needed to under- 
stand whether or not they are common, significant per- 
turbers of lens models. 

If many lenses do lie in groups, a comprehensive study 
of lens galaxy environments will provide a unique sample 
of galaxy groups with redshifts out to z ~ 1, the range 
of lens galaxy redshifts. The group environment is an 
important laboratory for studying galaxy evolution, as 
groups are common environments for galaxies and are 
relatively simple systems in which the range of factors 
impacting galaxy evolution (primarily galaxy-galaxy in- 



teractions) is much narro wer than in the complex, hot, 
dense cluster environment l)Zabludoff fc Mulchaevll99 8). 
There is a small but growing n umber of su rveys for 
groups at intermed i ate redshifts ijCarlberg et al. 200 5 
iWilman et al. "20051 iGerke et all '2005V These searches 
are challenging, because poor groups are difficult to find 
given their low projected surface densities and faint X- 
ray luminosities. The large, homogeneous sample of 
groups obtained from a lens survey would, when com- 
pared with local samples, permit direct observation of 
group evolution. 

One measure of group evolution may be the prop- 
erties of the brightest group galaxy (BGG). In 
nearby, high velocity dispersion (cr > 300 km s^^) 
group s, the BGG is typically a giant, central ellip- 
tical (|Zab]udoff fc Mulcha,e;il1998t Ivan den Bosch et all 
120051) . In contrast, low-cr groups at z ?» rarely have 
a dominant, bright member at their centers and tend 
to be elliptical-poor. These results suggest that higher- 
(7 groups are generally more dynamically evolved than 
their low-cr counterparts. It is unknown whether this di- 
chotomy persists out to higher redshifts or whether group 
evolution is sufficiently rapid to make distant groups with 
central giant ellipticals rare. A sample of poor groups at 
higher redshift may permit this question to be addressed 
for the first time. 

Because lensing select s on the basis of mass (e.g., 
iFukugita fc Turneil [T99l|) . one might first expect that 
lens galaxies in groups would also be the BGGs. How- 
ever, the sheer number of dwarf galaxies in a rich envi- 
ronme nt can offset the larger cros s section of the massive 
BGG l(Keeton. Christlein. fc Z"ab ludofi 2000), leading to 
the prediction that fewer than ~50% of lens galaxies are 
BGGs (Oguri 20061). Lens theory further suggests that 
lens galaxies span a large range of luminosities, with 
a typ ical luminosity around L* (e.g., iKochanek et alJ 
l200(l . The data are now in hand to determine directly 
the distribution of lens galaxy luminosities and the frac- 
tion of lens galaxies that are also BGGs. 

In order to address these issues, we have undertaken a 
photometric survey of 80 strong gravitational lenses, 
nearly the complete sample of lens systems known as of 
early 2004, with a significant fraction of these (~ 1/3) 
targeted for follow-up spectroscopy. The goals of our 
survey are to characterize the local environment of each 
lens galaxy out to large radii (~ 15 arcminutes), to iden- 
tify any line-of-sight structures toward each lens system, 
to determine the impact of the local environment and the 
interloping systems on the lens model, and to compile a 
large catalog of galaxy groups at intermediate redshifts 
useful for studying group evolution from z « 1 to the 
present. 

In IKeeton fc Zabludo3 l|2004(l . the theoretical under- 
pinnings of the biases introduced by over-dense environ- 
ments at the lens redshift are discussed in detail. In 
this paper we present the first results of our photometric 
survey, focusing on the fields of twelve lenses. Follow-up 
spectroscopy of eight of these lenses and detailed analyses 
of the shear and convergence introduced by both the im- 
mediate lens environment and any interloping structures 
are presented in MO©. In §2 of this paper we describe the 
lens sample, observations, and data reduction techniques. 
In §3 we discuss the red-sequence finding algorithm we 
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use to detect candidate groups in our photometric cat- 
alog, including simulations of its effectiveness and com- 
parisons to the spectroscopic data. In §4 we employ our 
sample of detected groups to address the issues raised 
above. Throughout this paper we assume a standard 
cosmology of Hq — 70 km s~-^ Mpc = 0.3, and 

f^A = 0.7. Unless explicitly stated otherwise, all physical 
quantities from referenced papers have been adapted to 
this cosmology. 

2. OBSERVATIONS 
2.1. Sample Selection 

We select our sample of twelve lens systems from the 
sample of all known lens systems on the basis of redshift 
(-^lens < 0.5) and foreground reddening {E{B — V) < 0.2). 
We also select roughly equal numbers of quadruple- and 
double-image lenses, as well as a handful of "other" im- 
age types (two objects with Einstein rings and one system 
with two lens planes). Imaging and analysis of the re- 
maining lens systems will be presented in a future paper. 
The current lens sample is given in Tabled 

Of these twelve lens systems, three have previously 
published spectroscopically confirmed groups at the 
lens redshift: MG J075H-2716 (hereafter MG0751; 
iTonrv fc Kochanekl IT9 99). PG 1115+080 (hereafter 
PG1115: lKundic et al."'1997a). and CLASS B1422+231 
(hereafter B1422; fKundic et al. 1997a). B0712+472 
(hereafter B0712) has a spectroscopic foreground group 
at z ~ 0.29 but no known group surrounding the 
lens (jFassnacht &: LubinI 1200^ . Two other lens sys- 
tems have either an over-density of galaxies projected 
around the len s (MG J1654+1346, hereafter MG1654; 
iLangston et al.l Il989() or a group identified by pho- 
tometric redshifts (HE 2149-2745 hereafter HE2149 ; 
iLopez. Wucknitz. fc Wisotzkil Il998t iFaure et alJl2004D . 
Two more lens systems are best characterized by mod- 
els including external shear, suggesting the presence of 
an additional mass compon ent: BRI 0952-0115 (here- 
after BR,I0952: lL'ehar et al.li200 0) and IR XS J11 3155.4- 
123155 (hereafter RXJ113i7E1usc et alj |2003|) . For 
the remaining four lenses, FBQS J0951-I-2635 (here- 
after FBQS0951), CLASS B1600-t-434 (hereafter B1600), 
PMN J2004-1349 (hereafter PMN2004), and CLASS 
B2114-I-022 (hereafter B2114), no evidence of a complex 
environment has been found previously. 

It is not yet clear how representative these twelve sys- 
tems are of all lens environments. This sample includes 
all eight lenses analyzed spectroscopically in M06, four 
of which were targeted there because they had known 
groups or good evidence for external shear. We use no 
information about environment to select the other four 
lenses in the present sample. 

2.2. Data collection 

We obtained observations over several observing runs 
between 2002 and 2004 with the Mosaic- 1 imager on the 
KPNO Mayall 4-m telescope and the Mosaic II imager 
on the CTIO Blanco 4-m telescope. The observing log 
is given in Table [3 The wide field of view of the Mosaic 
cameras (^ 36' x 36') covers the full virial radius, Tvir ~ 
500/i^^kpc ~ 3', of any groups that may be present 
around the lens galaxies, in addition to all other massiv e 
structures that may affect the lens potential (see lM06j) . 



We imaged each field with the "nearly Mould" / filter, 
and with either Harris V or Harris R, with the choice of 
bluer filter based upon the location of the 4000A break 
at the lens galaxy redshift. We changed from ^ to i? 
for zions > 0.35, approximately the redshift at which the 
4Q0QA break enters the R filter. ^ 

In order to fill the gaps between the CCDs in the Mo- 
saic cameras, we dithered the exposures. During the 2002 
runs, no specific dither pattern was used; in subsequent 
runs we used a five-point dither pattern. We based our 
exposure times in 2002 (~ 1800 sec in R and /) on those 
necessary to obtain S/N= 5 at / « 21.5. Subsequently, 
we found this signal-to-noise to be insufficient for star- 
galaxy separation for / > 21, and thus we increased the 
total exposure times to 2400-2700 sees. We split the 2002 
exposures into two 900 sec exposures in each band. This 
strategy proved less than ideal, as the signal-to-noise in 
the areas affected by the chip gaps was appreciably lower, 
some small gaps in areal coverage remained, and the I- 
band sky brightness was sufficiently high (^11000 ADU) 
to fill the pixel wells to ~30% capacity, appreciably re- 
ducing dynamic range. In subsequent observing runs, 
we split the total exposure time into 480 sec or 540 sec 
exposures. This improved strategy increased the signal- 
to-noise in the chip gap regions and kept the sky levels 
at reasonable values. 

2.3. Data Reduction 

We reduced the data using the MSCRED package of 
IRAF.^ Our data reduction closely followed the prescrip- 
tion used by the NO AO Deep Wide-Field Survey as out- 
lined in the online notes by Jannuzi, Claver, & Valdes.* 
In short, we combined each night's calibrations (bias 
frames, dome flats, and night-sky flats) and applied these 
to the science images. For KPNO data, we also created a 
map of the "pupil-ghost" image and subtracted this from 
each frame. For each night's /-band data, we created, 
scaled, and subtracted fringing maps from each frame. 
At this point, we used coordinates from the Guide Star 
Catalog 2 to refine the image world coordinate systems.^ 
The images were projected onto the tangent plane and 
then stacked using ia clipping. Unlike the Deep Wide- 
Field Survey reduction routine, we did not remove cos- 
mic rays in a separate step, as the a clipping removed 
the vast majority of hits. We also did not fit and remove 
sky gradients prior to the image stacking, as the resulting 
stacked images were generally cosmetically better with- 
out the sky removal. We set pixels in the final stacked 
image that had no good values (i.e., unobserved regions, 
saturated pixels, and bleed trails) to an arbitrary, high 
value for ease of flagging in the analysis. 

^ We have since modified this strategy such that our ongoing 
survey now uses V out to ziens = 0-5i ^■s the V filter gives better 
redshift resolution than the _R filter out to these larger redshifts. 
However, in the present paper, the R — I color proves sufficient to 
identify groups for 0.35 < ziens ^ 0-5, both l ocal to the lens and 
projected along the line-of-sight, as shown in il3.2l 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

^ http : //www . noao . edu/noao/noaodeep/Reduct lonOpt/f rames .html] 

^ The Guide Star Catalog was produced at the Space Telescope 
Science Institute under U.S. Government grant. These data are 
based on photographic data obtained using the Oschin Schmidt 
Telescope on Palomar Mountain and the UK Schmidt Telescope. 
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TABLE 1 

Selected lens systems 



Lens Name 






E{B~V) 


Obs. Nighf^ 




Exp. time (s) 


Calib. Nighf^ 




References'^ 










V 


R 


/ 


V 


R 


/ 


V 


R^ 


/ 




CLAbb B0712+472 


4 


0.414 


0.113 




K1,K4 


K3 




1800 


2700 








1,2 


MG J0751+2716 


R 


0.350 


0.034 




K2 


K3 




1800 


2700 








3,4 


FBQS J0951+2635 


2 


0.24^ 


0.022 


K2,K3 




K3,K4 


2700 




2700 








5,6 


BRI 0952-0115 


2 


0.41 


0.063 




C2 


C3 




1800 


1800 




C5 


C5 


6,7 


PG 1115+080 


4 


0.310 


0.041 


C3 




CI 


1800 




2400 


SDSS 






8,9,10 


IRXS J113155.4-123155 


4 


0.295 


0.035 


C7 




C7,C8 


2400 




2400 








11 


CLASS B1422+231 


4 


0.339 


0.048 


Kl 




K3 


1800 




2700 








9,10,12 


CLASS B1600+434 


2 


0.414 


0.013 




K2 


K4 




1800 


1800 






K8 


2,13 


MG J1654+1346 


R 


0.254 


0.061 


C2 




C1,C2 


2700 




4500 


K8 






14 


PMN J2004-1349 


2 


< 0.36 


0.202 




C2 


C3 




1800 


1800 




C9 


C9 


15,16 


CLASS B2114+022 


2+2? 


0.316, 0.588 


0.072 




C2 


C2 




1800 


1800 




C9 


C9 


17,18 


HE 2149-2745 


2 


0.495 


0.032 




CI 


CI 




1800 


1800 








19,20 



References. — fl) [ Jackson et alj|19 9a (21 IFassnacht fc Cohenlll 99a (3) iLehar et alj|1997t (41 ITonrv fc Kochaneklll99a (sySchechter irtalJ ll99a 
(6) 'Koch£mcForn?'2()0(in?^ 

[Patnaik ct aL 1992, (13) Jackson ct al. 1995, (14) Lansston ot al. 1989, (151 IWinn et alJl20nil (leT Twinn et alll20na (171 IKing fc BrownellTgM (18) 
"gusto ct aL..2001, (191 .Wisotzki et al 1990. (20l .Burud ct al..i2002. 

^Number of images; R indicates an Einstein Ring 
''From literature or CASTLES 
•^Scc Tabic |2l 
^Where necessary 

^References for lens discovery and lens galaxy rcdshift 
^Photometric Redshift 



TABLE 2 
Observing Log 



Night 


UT Date 


Telescope 


Instrument 


Seeing (") 


Comments 


Kl 


2002 March 14 


KPNO Mayall 4m 


Mosaic- 1 


1.3 




K2 


2002 March 15 




1.6 




K3 


2002 March 16 






1.2 




K4 


2002 March 17 






1.3 


cirrus; non- photometric 


K8 


2004 June 15 






1.0 




CI 


2002 May 10 


CTIO Blanco 4m 


Mosaic II 


0.8 - 1 




C2 


2002 May 11 






1.0 


non-photometric 


C3 


2002 May 12 






1.0-1.8 


non-photometric, volatile seeing 


C5 


2003 December 14 






0.8-1.0 




C7 


2003 December 18 






0.8-1.0 




C8 


2003 December 20 






0.8-1.0 




C9 


2004 June 18 






0.9-1.2 





Photometric measurements — We obtained the 
lens field photometry using SExtractor v. 2. 3. 2 
l|Bertin k ArnoutsI '1996'). Using the image in the 
best-seeing band (usually /) as a reference, we shifted 
and scaled the stacked images for a given field with the 
IRAF geotran task. We convolved the best-seeing image 
with a Gaussian kernel to degrade the seeing to that of 
the other band. 

Total magnitudes are iKronI l)1980|) magnitudes given 
by the SExtractor MAG_AUTO measured from the un- 
convolved /-band image. We measured galaxy colors via 
aperture magnitudes in the matched-seeing images. We 
selected the aperture size to be a 6.1 kpc physical ra- 
dius, equivalent to 5 pixels (= 1'.'3) for a z = 0.5 galaxy. 
This aperture corresponds to the full width-half max- 
imum (FWHM) of the point-spread function (PSF) in 
the worst-seeing lens field in this sample. 

Star-galaxy separation — We tested three methods for 
star-galaxy separation: (1) the stellarity index output by 
SExtractor, which compares observed profiles with an in- 



put FWHM value using a neural network algorithm, (2) 
measurements output by DAOphot (Stetson 1987), 
which are based on the goodness-of-fit of a model PSF, 
and (3) the difference in aperture magnitudes between 
two concentric apertures, with radii approximately equal 
to the FWHM and twice the FWHM. The last of these 
methods is motivated by the fact that point sources have 
a universal PSF and thus should have a constant dif- 
ference in aperture magnitudes, while extended objects 
(galaxies) should not. We inserted artificial stars and 
galaxies into the PG1115-h080 field using the IRAF art- 
data routine, recovered these objects with SExtractor, 
and performed star-galaxy separation using each of the 
three methods. 

The DAOphot analysis performed worst at star-galaxy 
separation due to the known high-frequency variations in 
the Mosaic PSF.^ The other two methods worked com- 
parably well over a wide range of parameters, but the 

^ See ||http://www.physics.nau.edu/~pmassey/Mosaicphot.htmH 
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Fig. 1. — Difference in aperture magnitudes for two concentric 
apertures as a function of magnitude in the field of B1422. Crosses 
indicate objects that are flagged as saturated or have otherwise 
bad photometry; filled circles are objects identified as stars; open 
symbols are galaxies; and filled triangles are objects that we then 
classified by eye. A visual inspection of the saturated objects finds 
only one bright galaxy near the edge of the field; the other objects 
are either saturated or bad measurements (e.g. near the image 
edge, under a bleed trail, a defect, etc.). The star-galaxy separation 
routine works well for 16.5 < / < 21. 



aperture method was able to go slightly fainter. 

The SExtractor stellarity index is a continuous pa- 
rameter running from 1 for point sources to for ob- 
viously extended sources, with intermediate values indi- 
cating various degrees of ambiguity. This parameter has 
two main drawbacks. First, it is quite sensitive to the 
input value of the FWHM for faint objects. Second, the 
value of the stellarity index corresponding to the dividing 
line between stars and galaxies is ill-defined. The divi- 
sion between stars and galaxies is subjective, and many 
different cuts exist in the literature. 

The concentric-aperture method is illustrated in Fig- 
ure n This method was robust over a large magni- 
tude range, though it failed at bright magnitudes due 
to an ambiguity between bright extended objects and 
saturated stars. Although SExtractor output saturation 
flags, some saturated stars failed to be flagged. There- 
fore, we classified by eye all non-flagged objects brighter 
than the magnitude at which saturation set in. For most 
fields, this involved fewer than a dozen objects. Due to its 
superior ability to identify unresolved objects, we there- 
fore used the concentric-aperture method for star-galaxy 
separation. 

2.4. Calibrations 

We calibrated our data obtained on photometric nights 
to the Kron-Cousins filt er set using M osaic images of 
Landolt standard fields l)Landomil992|) . We performed 
multiple-aperture photometry of photometric standards 
using DAOphot. Aperture corrections to infinity were 
calculated using the profile fitting routine DAOgrow. We 
assumed color transformations of the form: 

v = V + ao + ai{V-I)+a2iX~1.25) (1) 



r^R + bo + bi{R-I) + b2{X -1.25) (2) 
z = / + co + ci(l/-/) + C2(X-1.25), (3) 

where v, r, and i are the instrumental magnitudes and 
X is the airmass. In a few cases, depending on the bands 
observed during a night, the color term used in Eqn. ^ 
was B — V, and the color term used in Eqn. 13 was R — 
I; these are noted in Table 01 We assumed extinction 
values of 02 — 0.15, 62 — 0.10, and C2 = 0. 07 magnitudes 
per u nit airmass for the Kitt Peak data ijMassev et alJ 
I2002D . as we lacked sufficient standard star observations 
to determine these values. We were able to calculate 
extinction coefficients for each night of CTIO data. We 
determined color terms ai, bi and ci for each telescope 
run and made no attempt to account for known, small 
(millimag) chip-to-chip variations in color terms. The 
zero points oq, 60 and cq were determined individually 
for each night. The adopted photometric coefficients are 
given in Table |2| 

We applied the zero points and extinction coefficients 
to calibrate the galaxy photometry and colors but did not 
apply the color terms. The color terms derived above are 
for stellar spectral energy distributions (SEDs), which 
are roughly blackbody SEDs spanning a range of tem- 
peratures. The application of these color terms to colors 
of redshifted stellar systems was not appropriate, and at 
present we do not have calibrated spectrophotometry of 
any of our targets to determine empirical color terms for 
these galaxies. If we assume that these color terms are 
similar in size to those in Table 13 then the systematic 
errors we introduced into the colors would be on the or- 
der of 0.1 to 0.2 mags, similar to the size of the random 
errors ( W2.b\i . For galaxies with similar spectral energy 
distributions at a similar redshift, such as those in the 
red sequence (JJOJ, the systematic offset is the same 
for each galaxy and will not affect the detection of the 
red sequence. Therefore, we simply acknowledge that no 
color term correction has been applied to the photom- 
etry, likely introducing some small systematics into the 
derived colors (< 0.1 mags) and into the photometry- 
based group redshifts. 

In the cases where fields were imaged on non- 
photometric nights, we obtained short calibration expo- 
sures of the fields on later photometric nights. We de- 
rived local standard stars in each field using aperture 
photometry, including appropriate aperture corrections. 
From these local standards we determined zero points, 
which include the atmospheric extinction implicitly, for 
each field. The fields calibrated in this manner, and the 
photometric nights used for calibration, are indicated in 
Table n 

Special calibration steps were required for the V- 
band imaging of PG1115, for which we failed to obtain 
an image during a photometric night. Our attempts 
to calibrate the field with archival HST and archival 
CFHT data failed due to a lack of common unsaturated 
stars between the archival data and our field. Luck- 
ily, the PG1115 field was included in the Third Data 
Release of the Sloan Digital Sky Survey (SDSS DR3, 
Abazaiian et al.. .2005.) .^ We retrieved PSF magnitudes 
for all stellar objects within a 15' radius of PG1115 from 

^ Funding for the Sloan Digital Sky Survey (SDSS) has been 
provided by the Alfred P. Sloan Foundation, the Participating In- 
stitutions, the National Aeronautics and Space Administration, 
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TABLE 3 
Adopted Photometric Coefficients 



Night'' 


Band 






Fit Coefficients'^ 






N 


Kl 


V 


ao = 


-25.163 ±0.002 




n 034 + 002 




a2 = 0.15° 


4 




R 


bo = 


-25.369 ±0.005 




—0.052 ± 0.006 




b2 = 0.10° 


4 




I 


CO = 


-24.774 ± 0.010 




—0.013 ± 0.004 




C2 = 0.07° 


3 


K2 


V 


ao = 


-25.101 ± 0.007 


(3-1 — 


n 034 + n no2 




a2 = 0.15° 


15 




R 


bo = 


-25.390 ±0.006 


h-\ — 


— 052 + 006 




b2 = 0.10° 


10 




I 


CO = 


-24.777 ±0.003 


Cl — 


— 01 3 + 004 




C2 = 0.07° 


11 


K3 


V 


ao = 


-25.142 ± 0.004 




0.034 ± 0.002 




a2 = 0.15° 


4 




R 


bo = 


-25.343 ± 0.003 


'-'1 


—0.052 ± 0.006 




b2 = 0.10° 


6 




I 


CO = 


-24.777 ± 0.003 


Cl — 


— 01 3 + 004 




62 = 0.07° 


11 


K8 


V 


ao = 


-25.077 ± 0.023 


fll = 


0.001 ± 0.023 




a2 = 0.15° 


18 




R 


bo = 


-25.303 ± 0.021 


bi = 


-0.081 ± 0.040 




62 = 0.10° 


13 




I 


CO = 


-24.732 ± 0.017 


Cl = 


-0.031 ±0.016 




C2 = 0.07° 


14 


CI 


R 


bo = 


-25.667 ± 0.011 


bi = 


0.013 ±0.027 


b2 


= 0.031 ± 0.045 


8 




I 


CO = 


—24.912 ± 0.005 


Cl = 


0.027 ± 0.012=^ 


C2 


= 0.000 ± 0.022 


10 


C5 


V 


ao = 


-25.330 ±0.003 


ai = 


0.054 ± 0.005 


a2 


= 0.145 ±0.006 


72 




R 


bo = 


-25.517 ±0.003 


bi = 


-0.022 ±0.006 


b2 


= 0.097 ±0.006 


77 




I 


CO = 


-24.860 ± 0.004 


Cl = 


0.028 ± 0.004 


C2 


= 0.076 ± 0.007 


73 


C7 


V 


ao = 


-25.307 ±0.004 


ai = 


0.036 ± 0.003° 


a2 


= 0.068 ± 0.036 


58 




R 


bo = 


-25.489 ± 0.004 


bi = 


-0.010 ±0.005 


b2 


= 0.087 ± 0.049 


26 




I 


CO = 


-24.829 ± 0.006 


Cl = 


0.027 ±0.004 


C2 


= 0.065 ± 0.072 


46 


C8 


V 


ao = 


-25.341 ±0.008 


ai = 


0.043 ± 0.007° 


12 


= 0.135 ±0.024 


39 




R 


bo = 


-25.521 ±0.006 


bi = 


-0.029 ±0.010 


b2 


= 0.092 ±0.015 


32 




I 


CO = 


-24.872 ± 0.006 


Cl = 


0.028 ± 0.005 


C2 


= 0.046 ±0.017 


53 


C9 


R 


bo = 


-25.610 ±0.006 


bi = 


-0.015 ±0.013 




b2 = 0.090° 


45 




I 


CO = 


-24.963 ± 0.007 


Cl = 


0.070 ± 0.015'^ 




C2 = 0.050° 


23 



^Photometric nights only 
''See Eqns.IHtoEl 

° Adopted from published mean extinction 

Coefficient for ij - / 
''Coefficient for B-V 



the SDSS DR3 database. We transformed these magni- 
tudes to UBVRnln using the observed transformations 
given m iSmith et a l. ( 2002); we then used these magni- 
tudes as local standards in our V and / images. The 
resulting zero-points were oq = 25.390 ± 0.018 for the 
y-band and cq = 24.947 ± 0.052 for the /-band. The 
/-band zero point was consistent with that derived from 
the Landolt standards, giving us confidence that the V- 
band calibration was accurate. 

2.5. Completeness and Photometric Accuracy 

We determined the completeness and photometric ac- 
curacy of the observations via artificial galaxy tests. Ar- 
tificial galaxies were placed at random locations in the 
field, with additional artificial galaxies placed at random 
within Tvir = 500ft.^^kpc of the lens galaxy. The input 
galaxy colors were chosen to match those expected for an 
early-type galaxy at the lens redshift (see Appendix A). 
The input magnitudes were distributed randomly and 
equally between /* — 2 and / = 22.5, where /* is the 
observed apparent / magnitude of an galaxy located 
at zions (see Appendix A). For the lenses in this paper, 
17.40 < /, < 19.26. This range includes the bulk of the 
galaxies in our photometric catalog. 

Galaxy luminosity profiles were de Vaucouleurs {r^/^) 
profiles convolved with a Gaussian kernel whose FWHM 
matched that measured in the image. We calculated an 

the National Science Foundation, the U.S. Department of Energy, 
the Japanese Monbukagakusho, and the Max Planck Society. The 
SDSS Web site is http://www.sdss.org/ 



effective radius from the Kormendy relation 

(^)e =a + /31ogi?e, (4) 

where Re is the effective radius in kpc and (/i)e is the 
average surface brightness int erior to Rp in uni ts of mag- 
nitudes per square arcsecond <)Kormendvll977(l . By defi- 
nition, the total light within R^ is half of the total galaxy 
luminosity, so 

= 2.5 log(27r) + 5 log(re) +mT, (5) 

where r^, is the effective radius in units of arcseconds and 
mT is the total apparent magnitude of the galaxy. Com- 
bining these two e quations, setting /? = 3 (the approx- 
imate value from iLa Barbera et al]l2003() . and solving 
for logi?e, we get 

logi?, = I - 1.251og(27r) - ^log^kpc - ^ , (6) 

where 6'kpc is the angle (in arcsecond s) subtended by 1 
kpc at the lens galaxy redshift. For of. lLa Barbera et alJ 
ll2003|) give a rest-frame F-band value of 

av = 19.05 + 7.7 log(l + Zjens) + Ay , (7) 

where Ay is the F-band extinction. We converted this 
value to the observed /-band with a /^-correction com- 
puted with the synphot tools in IRAF's stsdas p ackage, 
using the E/SO spectrum from iColeman et al.l (^80). 
We also applied Galactic reddening corrections interpo- 
lated from the Schlegcl et al. ( 1998) maps using the red- 
dening curve of iRieke fc Lebofskvl ((19850 ■ 

We inserted the resulting artificial galaxies into each 
lens field and recovered them using SExtractor. Each run 
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included 365 artificial galaxies: 350 randomly scattered 
throughout the field and an additional 15 within r^ir of 
the lens. We considered an artificial galaxy recovered if 
its centroid was within 2.5 pixels and its / magnitude and 
color were within 1.0 mag of the input values. Fifty runs 
were performed in each lens field to build up statistics. 

We plot the results from the artificial galaxy tests in 
a typical field (B1422) in Figure |2] Two systematic is- 
sues are readily observed. The recovered galaxy mag- 
nitudes consistently underestimated the flux, the effect 
being greatest for the brightest galaxies. This underes- 
timate is due to the low surface brightness outer regions 
of an elliptical galaxy not being detected by SExtractor. 
The effect is strongest at bright luminosities because in- 
trinsically bright ellipticals have the most diffuse surface 
brightness profiles. The net effect is to amplify the num- 
ber of galaxies around /, and suppress the number of 
galaxies brighter than ~ — 1. For galaxies in these 
magnitude ranges, another method of measuring magni- 
tudes (such as profile fitting) is necessary to accurately 
determine the total luminosity. Nevertheless, we did not 
attempt to correct the photometric catalog for these sys- 
tematic magnitude offsets for two reasons. First, our 
method of measuring magnitudes is in common use; cor- 
recting to total magnitudes would hinder comparisons 
with other samples. Second, the exact offset depends 
on the input galaxy luminosity profile. Our tests were 
limited to bulge-dominated galaxies, so we could not es- 
timate the corrections for disk-dominated galaxies. In 
future work we will perform detailed profile fits to the 
galaxies near the lens systems to address this systematic 
effect carefully. 

For the brightest bins in the figure, the completeness 
is low because the central regions of bright galaxies are 
flagged as saturated. In order to determine if this sat- 
uration affects our results, we visually inspect objects 
flagged as saturated in a subset of the fields. While a 
handful of bright galaxies in each field are flagged as 
saturated, in no case are the flagged galaxies potential 
group members; these galaxies are either more than sev- 
eral arcminutes from the lens or obviously foreground, 
with magnitudes significantly brighter and physical ex- 
tents significantly larger than any other galaxy near the 
lens. 

Given that galaxies in groups tend to be clustered, 
we test whether crowding (i.e., the blending of overlap- 
ping galaxies) could be a significant source of incomplete- 
ness in our data. Overall, we find that the completeness 
within Tvir of the lens is similar to that in the remainder 
of the field. In some fields, such as B1422 (Figure[2Il, we 
find the completeness within rvir of the lens is slightly 
higher than in the field as a whole. In other fields the 
completeness is slightly lower. These small completeness 
variations are most likely due to statistical fluctuations 
resulting from the smaller number of artificial galaxies 
used to calculate the completeness in the relatively small 
area surrounding the lens as compared to the much larger 
area of the complete Mosaic field. We therefore conclude 
that crowding is not an issue in our data. 

Given the uncertainties in the model galaxies, the (ig- 
nored) scatter in the Kormendy relation, and the fact 
that the input artificial galaxies are solely early-type 
galaxies at a fixed redshift, the completeness is not nec- 
essarily valid outside the red sequence and is subject to 



uncertain systematics. However, the observed field-to- 
field stability of the completeness, despite differing lens 
redshifts, seeing, and field depth, suggests that our com- 
pleteness calculations are relatively robust. 

3. IDENTIFYING GROUPS AND CANDIDATE GROUP 
MEMBERS 

Due to the morphology-density relation, the early-type 
fraction is higher in overdense regions like groups and 
clusters of galaxies than elsewhere. In such cases, one 
observes a "red sequence" (RS) of elliptical and SO galax- 
ies in the color-magnitude diagram (CMD) of galaxies 
in the region. The RS forms the basis for our success- 
ful two-band photometric detection of galaxy groups and 
clusters. 

Our RS detecti on algorithm is a modi fied version of 
that presented bv iGladders fc Yed l)2000D . Their algo- 
rithm searches color and astrometric space for localized 
enhancements of early-type galaxies at similar redshifts 
and has proven successful in locating group-mass struc- 
tures over a wide range in redshift using data from op- 
tical imaging surveys. As our survey is targeted toward 
the immediate environment of lenses, we do not need 
to include the spatial filtering in our selection. The 
other major diffe rence between our algorithm and that of 
Gladders & Yei l|200(i>) is that we include human inter- 
action in the selection process. While this introduces 
a hard-to-quantify bias into the resulting group sam- 
ple, it enables us to detect less-massive groups (cr < 
300 km s~ ) . As shown in ii3.3l the resulting group cata- 
log is nearly complete with few false detections, suggest- 
ing that any bias introduced by the human interaction is 
small. 

3.1. Group Detection 

Our basic RS detection algorithm is as follows. First 
we search within rvir of the lens galaxy for an excess 
of galaxies with similar colors above the measured back- 
ground color distribution. We fit an RS to these galaxies, 
and calculate the projected spatial centroid of the can- 
didate group. Using the RS fit and centroid, we repeat 
this search to ensure that the group is not spurious and 
to refine the RS parameters and centroid measurement. 
We continue the iterative process until the changes in 
the centroid and the RS parameters are within the cal- 
culated errors. The first time we apply this algorithm to 
a field, we focus on galaxy colors near that of an early- 
type galaxy at the lens redshift, zjons- We then search 
for line-of-sight groups by repeating the algorithm until 
we have investigated all peaks in the initial color distri- 
bution. 

We assume an RS of the form 

iV-I)o^iV-I).,-l3viIo-h), (8) 

where /3y is the slope of the RS, and {V — /)* is the 
unreddened observed V — I color of an galaxy at the 
lens redshift. For fields with i?-band imaging, we use this 
same relation with V replaced by R. 

We allow both the RS slope (3v and the group red- 
shift, which determines /* and {V — /)», to vary freely. 
We calculate the initial {V — I)^. and J* from the galaxy 
models described in Appendix A and the known zions- 
The RS slope has bee n found to have values around d(U— 
V)/dMv = 0.08 re.g.. lBower et al.ll992tlBeU et al.l2004 
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Fig. 2. — Results of artificial galaxy tests in the field of CLASS B1422+231 (hereafter B1422). The top panel shows the number of 
artificial galaxies (over all 50 runs) as a function of input total magnitude for all galaxies (heavy solid line) and all recovered galaxies 
(dotted line) , as well as the number of recovered galaxies as a function of measured magnitude (dashed line) . The second panel shows the 
completeness as a function of input total magnitude for the entire image (heavy solid line) and the region within a 500h~^ kpc projected 
radius of the lens (dashed line). The third panel shows the mean difference between the total magnitude and the measured magnitude as 
a function of total magnitude; error bars give the la deviation about the mean. The bottom panel shows the mean difference between the 
input and measured aperture colors as a function of total magnitude; error bars are the Icr deviation about the mean. 



iMcIntosh et al.l2005D . Because the color-magnitude rela- 
tion is difficult to model properly, and because our band- 
passes do not correspond exactly to rest-frame U and V, 
it does not necessarily follow that this slope can be di- 
rectly compared with the V — I ot R — I slopes of the 
groups in this paper. However, we show in ii4.1.1l that 
the actual fit slopes from our sample are indistinguish- 
able from this value. We therefore initially set /3y = 0.08. 
Although this could bias the final fit slopes, we show in 
21 that we can recover the correct RS slope even if this 
initial assumed value is incorrect. 

For the first iteration of our RS search, we limit the 
search to within rvir (= 500h~^ kpc) of the lens galaxy. 
This radius is large enough (~ 2'-3') to include most of 
the line-of-sight structures that could impact lens models 



(see IM06I) . It is possible that a massive galaxy cluster 
could lie outside this radius and still alter the lens po- 
tential, so we estimate the fraction of lenses affected by 
such a cluster as follows. We consider a lens model per- 
turbation signi fican t if the convergence n > 0.05. From 
Eqn. A20 in iMOfiL we find that the maximum impact 
parameter 6niax that a cluster can have and significantly 
impact the lens model is 6max = 4.8 x 10~^a^, where b^^x 
is in arcminutes and a is the cluster velocity dispersion 
in km s~^. This conservatively assumes that the cluster 
is at the lens redshift; for clusters at other redshifts, 6max 
will be smaller. For a galaxy cluster with cr ~ 750 km 
s~^, the median velocity dispersion for rich Abell clus- 
ters {^Zabludoff. Huchra. & GcUcr 1990), 6niax « 3' — 
close to the radius within which we search anyway. The 
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surface density of gal axy clusters is ~ 11 per sq. degree 
ijGonzalez et aUEoOlL 0.3 < ^duster < 0.9, correcting for 
their 30% false detection rate). On average, this results 
in 0.07 galaxy clusters within fomax- We consider this to 
be an upper limit, given that we are already searching for 
groups and clusters out to a significant fraction of 6max 
and that 6max is smaller if there is any redshift difference 
between the cluster and the lens. We therefore conclude 
it is unlikely that any significant perturbing structures 
lie outside our search radius. 

We correct galaxy magnitudes for interstellar extinc- 
tion and reddening usin g standard Galactic extinction 
(|2,ieke & Lebofskv"1983) with Ry = 3.1 and extinctions 
from Schlegel et al. (1998) obtained from the NASA Ex- 
tragalactic Database.® We do not attempt any correc- 
tions for differential Galactic extinction across the field; 
any variations within rvir of the lens are likely to be neg- 
ligible at the Galactic latitudes of these lenses. As is 
common practice, extinction-corrected magnitudes and 
colors are denoted by the subscript "0." 

To account for the slope of the RS, we define an "ef- 
fective color" for each galaxy: 

iV-lU = iV-I)o + l3v{Io-h)- (9) 

This represents the color each galaxy would have if 
it resided on the RS and had an apparent magnitude 
/ = This definition is useful because all RS galaxies 
in a group should have the same effective color. Hence, 
we can identify RSs as peaks in a histogram of effective 
color, providing an efhcient way of searching for groups of 
early-type galaxies. The search for peaks in the effective 
color histograms is also adept at identifying groups with 
high late-type fractions, as blue galaxies in a group will 
all have roughly similar colors, although in the absence of 
morphological or spectroscopic data to classify the galax- 
ies as late-types, the derived redshift of the group will be 
incorrect (see i)3.2|l . 

We compute effective color histograms for galaxies 
within Tvir of the lens galaxy, as well as for all galaxies 
in the image (normalized to the area within rvir). The 
number density of field galaxies as a function of appar- 
ent magnitude has a ste eper slope tha n typical gr oup lu- 
mino sity functions (e.g. JLin et al.llT99 9: Gladdcrs fc Yed 
12000(1 . so it is possible to choose a magnitude limit that 
maximizes the contrast between the RS and the field 
galaxy distribution. This magnitude limit will vary from 
field to field depending on the richness of the group 
and statistical fluctuations in the number of background 
galaxies. Our empirical testing finds that, in most cases, 
the best contrast between the group and the field is 
achieved for a limiting magnitude of I ^ I^, + 2.5. We 
therefore consider only galaxies brighter than -|- 2.5 
throughout this algorithm. 

We visually examine the effective histograms for an ex- 
cess of galaxies of a similar effective colors relative to the 
normalized field background, an indication of a possible 
group. We select all galaxies within a distance A{V—I)cS 
of the peak in the effective color histogram as candi- 
date RS galaxies. We take A{V — I)cS to be the larger 

® This research has made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion Labo- 
ratory, California Institute of Technology, under contract with the 
National Aeronautics and Space Administration. 



of the photometric color error and 0. 05 mags, the typi- 
cal RS scatter in the literature fe.er.. iBower et alJll992t 
iMcIntosh et al.ll2005|) . The photometric color errors are 
typically larger than the intrinsic scatter for / > 21. 

We use linear regression to fit a linear RS to the se- 
lected galaxies brighter than /* -I- 2.5, thereby refining 
the RS slope (3v and intercept. From this fit we deter- 
mine the new values of {V — /)*, /*, and the RS redshift 
zrs based on the early-type galaxy photometric models 
(Appendix A). As discussed in Appendix A, we note that 
changes in the assumed normalization of results in sys- 
tematic changes in zrs , though for the range of published 
values of L*, this bias is Az^g < 0.02. We also calculate 
the raw and luminosity-weighted group centroids of the 
RS galaxies. 

At this point, we can either accept or reject the new 
parameter values and then choose whether another iter- 
ation should be performed. Changes in redshift cause 
the projected angular size of rvir to be recomputed, and 
changes in the centroids adjust the area of sky in which 
galaxies are examined. The iterative process continues 
until the change in the RS redshift between iterations is 
less than Azrs = 0.02 and the shift in the centroid is less 
than 1(7. If multiple peaks in the color distribution are 
observed, we repeat the entire process on each potential 
RS. For the purpose of this paper, we define "significant" 
candidate RSs as those with at least three selected galax- 
ies above the averaged background persisting for at least 
two iterations, criteria that produce the best results in 
tests with mock groups (see ^3.2|l . 

An example of this process applied to the lens B1422 
is shown in Figure |3| We see an RS of galaxies at 
V — I « 1.9, roughly the effective color of early- type 
galaxies at zjons and also the approximate color of the 
lens galaxy. At brighter magnitudes, the group lumi- 
nosity function clearly dominates that of the field. A 
map of the galaxies in this range (Fig. IHt) shows that 
the collection of bright galaxies, and therefore the group, 
is nearly centered on the lens galaxy and is elongated. 
This figure shows that our group detection algorithm 
is capable of detecting physical groups and constrain- 
ing the group redshift, richness, centroid, and shape. We 
note that some of the bright galaxies in the selected RS 
may not be group members but rather higher-redshift 
blue galaxies. In the future, we will also perform a mor- 
phological analysis, which should separate background 
disk-dominated galaxies from the early-type RS group 
members. For now, we make two tests of the algorithm, 
one using simulated groups of galaxies f ii3.2|l . the other 
using spectroscopy of eight of the twelve lens fields to 
verify RS membership ft l3.3(l . 

3.2. Testing the Algorithm: Simulated Data 

We test our RS search algorithm by inserting mock 
galaxy groups into each field's galaxy catalog. We con- 
struct mock photometric groups from empirical luminos- 
ity functions. The early-type fraction and number of 
bright group galaxies are readily modified to simulate 
galaxy groups of different richnesses. The precise lumi- 
nosity function does not significantly affect the success 
of the algorithm. 

We begin by creating a mock group galaxy catalog de- 
signed to cover a range of redshifts (z = 0.2, 0.4, 0.6), 
early-type fractions (/c — 0.0, 0.5, 1.0), RS slopes 
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Fig. 3. — Illustration of our group-finding algorithm for B1422. (a) Histograms of effective colors for galaxies brighter than 7» +2.5 
projected within rvir (solid) and for the entire field normalized to the area within r^ir (dashed). Vertical dotted lines bracket the region 
in (V — /)cff from which the final group member catalog was selected, (b, top panel) Color-magnitude diagram for all galaxies (crosses) 
projected within rvir. The star with error bars indicates the lens galaxy; concentric circles indicate the expected location of an L» early-type 
galaxy at the lens redshift. (b, middle panel) Color-magnitude diagram for galaxies projected within rvir (crosses) with galaxies selected 
by the red sequence finding algorithm shown as filled boxes. The solid line is the fit red sequence. The vertical dashed line in all right- 
hand panels indicates the faint magnitude limit of the final red-sequence fitting, (b, bottom panel) Residuals of the selected red sequence 
candidates about the best-fit line. Error bars indicate errors estimated from the artificial galaxy tests, (c) Map of red sequence-selected 
galaxies projected within rvir ^ ^ function of distance from the lens galaxy (open star). Filled circles indicate galaxies brighter than /« -f 2.5 
and are larger for brighter objects. This figure illustrates the ability of the red sequence finding algorithm to locate physical groups. 



{(3v =0.05, 0.10, 0.15), and richnesses. We define the 
mock group richness by the number of galaxies brighter 
than -|- 2, and our mock groups are designed to have 
richnesses of 4, 8, or 16. A richness of 4 corresponds to 
Local Group analogues. We derive cumulative luminosity 
functions fro m the empirical group l umino sity functions 
presented in lZabhidoff fc MulchaevI ()200nD to determine 
that richnesses of 8 and 16 correspond to velocity disper- 
sions of CTv ^ 300 km s^^ and CTv ~ 500kms~^, respec- 
tively, though there is significant scatter in this relation. 

We select the luminosity function for each group from 
the e mpirical l uminosity functions derived from GEMS 
data ijMiles et a l. 2004) . That study presents different lu- 
minosity functions for X-ray bright {L^ > 10*^ '' ergs~^) 
and X-ray faint /undetected groups. This luminosity 
bre ak corresponds to a velocity dispersion of ^ 200 km 
s-i l|Mulchaev fc ZabludollTool . so we use the GEMS 



empirical X-ray faint luminosity function for groups with 
richness 4 and the X-ray bright luminosity function for 
groups with richness 8 or 16. 

We assign each mock galaxy an / magnitude drawn 
randomly from a distribution defined by the appropri- 
ate luminosity function. We then label each galaxy an 
early- or late-type with a probability based on the input 
early-type fraction (/c). Late-type galaxies are further 
divided with equal probability into Sa, Sb, and Sc types. 
We assign colors based on the appropriate photometric 
model (Appendix A) for the assigned morphology, and, 
for early- types, we place the galaxy on the RS with a 
selected RS slope between /3y — 0.05 and 0.15. We add 
Gaussian scatter of ay-i = 0.05 mag, similar to the ob- 
served scatter in other RS studies, to each mock galaxy's 
color. We also add additional scatter to the / magnitude 
and color of each galaxy to simulate the photometric er- 
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rors based on the photometric accuracy determined from 
the artificial galaxy tests ( >i2.5(l . 

We now need to assign positional data to each galaxy. 
To simulate the foreground and background galaxies in 
the field, we assign the cataloged galaxies random x- and 
y-positions in the B1422 field (chosen to be a representa- 
tive lens field), thereby mixing any actual galaxy groups 
and clusters into the field galaxy population. We assume 
that all mock group galaxies lie within rvir of the lens, 
and we assign these galaxies positions drawn randomly 
from a 2-dimensional Gaussian density function centered 
on the lens position. We then run the group finding al- 
gorithm on the combined catalog. We consider a mock 
group to be recovered if we find a candidate RS with 
an effective color within 0.2 of the input group's color, 
a value similar to our input color errors. This criterion 
prevents the detection of spurious groups, i.e., those con- 
taining none of the input mock galaxies. We make an ex- 
ception to this criterion if the group has a high late-type 
fraction, in which case there often is a false RS created 
by having many blue galaxies with similar colors. In this 
circumstance, the group is considered recovered despite 
having an effective color >0.2 mag bluer than that of a 
true RS at the group redshift. In group catalogs based 
on our actual data, groups with high late-type fractions 
would be labeled line-of-sight groups. 

Our knowledge of the mock group color could poten- 
tially bias the group finding algorithm toward detection 
of the mock group. However, we also know the redshifts 
(and therefore expected RS colors) of the lenses in the 
present sample, so any bias for detecting groups at zions 
and the mock groups is the same. The small number of 
interloping groups identified from spec troscopy that are 
undetected photometrically (see argues that this 

bias is also small for line-of-sight group detection. The 
results of the simulations are given in Table El 

These tests indicate that we can detect rich groups 
via our detection algorithm out to the highest lens red- 
shifts in our sample, independent of /o. Groups with 
velocity dispersions of ~300 km s~'^ are detectable out 
to z « 0.4, although by z = 0.6 the detection rate is quite 
low (< 50%). Local Group analogs are detectable at the 
~50% level at low redshifts (z = 0.2) if /c is close to one 
or to zero, although the former case is unlikely to occur 
given that the early-type fractions of th e poorest groups 
are small (Zabludoff & MulchaevI 1 1998(1 . The detection 
rate of the Local Group analogues at z > 0.4 is consistent 
with the false-detection rate (see below). We are there- 
fore confident that we can reliably detect the majority 
of galaxy groups with velocity dispersions > 200-300 km 
s~^ out to z = 0.4. We can detect poorer groups out to 
z « 0.2. 

From the simulations we learn that, in the absence 
of morphological information, the fit RS colors are sys- 
tematically bluer than the input value. The reason is 
that some late-type galaxies (typically Sa's) are included 
in the sample when the RS is fit. For systems with /c 
near 1, the error is small. In the groups with low /c, we 
observe a broad peak in the effective color distribution 
due to the late-type galaxies, and the estimated group 
redshift is biased low (the magnitude of this bias is red- 
shift dependent). We could lessen this systematic effect 
by including morphological information in our galaxy se- 
lection; in future work we will use galaxy profile fitting 



software (e.g., G/MgD : Isiinard et al.l2'00^ to reject disk- 
dominated systems from the sample used to fit the RS. 
Even without this morphological information, we deter- 
mine group redshifts with reasonable accuracy (see ii3.3|l . 

We note that the individual fit RS slopes are very noisy 
and imprecise, although the ensemble average is rela- 
tively accurate. For groups with /c « 1, the individual 
slope measurements are much more accurate (Table QJ. 
This suggests that the accuracy of the slope determina- 
tions would be greatly improved by using morphological 
data to exclude late- type galaxies from the RS fit. 

Finally, we note that the measured scatter about the 
RS is not a reliable measure of the input scatter. For 
groups with J* -t- 2.5 > 20 (z > 0.25), the photomet- 
ric color errors dominate the RS scatter. By contrast, 
for groups with few RS members, the paucity of red se- 
quence galaxies results in an artificially low scatter. Only 
for richness 16 groups at low redshifts (z < 0.25) is the 
recovered scatter comparable to the input scatter. 

False positives — To determine the rate of false pos- 
itives, we run the group-finding algorithm on two lens 
fields, B1422 (for V - I) and MG0751 (for R-I). We 
assign each galaxy in the photometric catalog a new, ran- 
dom position in the field, thereby washing out any true 
groups and clusters. We then search for peaks in the ef- 
fective color histogram, both within 0.05 mags of the lens 
galaxy color (the typical scatter about the input color in 
the mock RSs), simulating false positive detections of 
lens groups, and at other colors, simulating false line-of- 
sight structures. Fifty such runs are performed in each 
field. 

For the V — I data, we record only two (i.e., 4%) false 
positive detections of lens galaxy groups and eight (16%) 
false positive line-of-sight groups. For the R — I data, 12 
fields (24%) have false positive lens galaxy groups, and 
six (12%) have false positive line-of-sight groups. The 
higher number of false positives in the R — I data arises 
mainly from the smaller color range spanned by the cat- 
alog galaxies (~2 mags, compared to ^^4 mags in V — I). 
Each color bin therefore contains more field galaxies, and 
the Poisson fiuctuations often exceed the group thresh- 
old. The rate of false positive lens groups in the MG0751 
field is also enhanced by the fact that the expected RS 
color at zicns [R — I — 0.91) lies near the peak of the field 
galaxy color distribution, so shot noise from the field pop- 
ulation is even larger. This latter problem vanishes for 
i?-band lens fields where (i? — /)* > 1, or ziens ^ 0.45, as 
the field density of galaxies at these colors is much lower. 

As a whole, false positives are not a serious concern for 
the lenses in our sample with V band imaging. For lenses 
with _R-band imaging, we need to account for the false 
positive rate, especially for lower-redshift lenses where 
the R — I colors are similar to that of the majority of 
background galaxies. As will be shown in M.l.ll these 
false positives have a minimal impact on our conclusions. 

3.3. Testing the Algorithm: Spectroscopic Confirmation 

The ultimate test of the group finding algorithm 
is whether the candidate groups can be confirmed 
as real, bound systems. We have obtained spectro- 
scopic redshifts of ~ 100 galaxies in each of eight lens 
fields (MG0751, BRI0952, PG1115, B1422, MG1654, 
PMN2004, B2114, and HE2149); the identification of 
bound systems and the analysis of their impact on lens- 



12 



Williams et al. 



TABLE 4 

Results from group-finding algorithm tests with mock groups 



/c 




Pv 


Mrs 


N It+2,det 




A{V-/)..flt 




Scatter 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




{°) 








1.00 


16 


0.05 


5 


13.3 ±3.2 


0.023 ±0.022 


-0, 


■ Oil ± 0.015 


0, 


,065 ± 0, 


,019 






0.10 


5 


14.1 ± 2.4 


0.005 ±0.014 


-0, 


.027 ± 0.019 


0, 


,080 ± 0, 


,015 






0.15 


5 


14.4 ±2.5 


-0.014 ±0.011 


-0, 


.022 ± 0.012 


0, 


,067 ± 0, 


,025 




8 


0.05 


4 


6.5 ± 1.0 


0.017 ±0.044 


0, 


.008 ± 0.027 


0, 


,079 ± 0, 


,016 






0.10 


3 


6.8 ±3.2 


-0.038 ± 0.046 


0, 


.001 ± 0.029 


0, 


,008 ± 0, 


,002 






0.15 


4 


6.0 ± 1.6 


0.009 ± 0.043 


-0, 


.007 ± 0.013 


0, 


,005 ± 0, 


,002 




4 


0.05 


1 


10.6 


0.028 




0.030 




0.007 








0.10 


3 


3.9 ± 1.1 


-0.026 ±0.032 


-0, 


.036 ± 0.015 


0, 


,097 ± 0, 


,012 






0.15 



















0.00 


16 




5 


10.7 ±4.8 


0.062 ±0.053 


-0, 


.184 ± 0.013 


0, 


,010 ± 0, 


,005 




8 




4 


5.2 ± 1.0 


0.047 ±0.034 


-0, 


.179 ± 0.020 


0, 


,013 ± 0, 


,006 




4 




3 


3.2 ± 1.0 


0.104 ±0.020 


-0, 


.186 ± 0.008 


0, 


,019 ±0, 


,004 


0.50 


16 


0.10 


5 


7.8 ± 2.9 


-0.001 ± 0.008 


-0, 


.047 ± 0.024 


0, 


,011±0, 


,004 




g 


0.10 


3 


5.4 ± 0.9 




-0, 


.052 ± 0.043 


0, 


,017 ± 0, 


,004 




4 


0.10 


1 


4.6 


0.012 




0.030 




0.011 




1.00 


16 


0.05 


5 


15.2 ±3.9 


0.026 ±0.023 


-0, 


.017 ± 0.016 


0, 


,107 ± 0, 


,017 






0.10 


5 


17.6 ±3.1 


0.000 ±0.010 


-0, 


■ Oil ± 0.016 


0, 


,087 ± 0, 


,027 






0.15 


5 


16.2 ±2.2 


-0.001 ±0.027 


-0, 


.023 ± 0.013 


0, 


,088 ± 0, 


,032 




8 


0.10 


5 


7.9 ± 2.2 


0.017 ±0.027 


-0, 


.017 ± 0.016 


0, 


118 ± 0, 


,026 




4 


0.10 


3 


5.2 ±3.1 


0.024 ± 0.054 


-0, 


.031 ± 0.054 


0, 


146 ± 0, 


,014 


0.50 


16 


0.10 


5 


8.4 ±2.3 


0.034 ± 0.043 


-0, 


.178 ± 0.219 


0, 


102 ± 0, 


,026 




8 


0.10 


4 


5.3 ± 1.8 


0.000 ± 0.067 


-0, 


,244 ± 0.249 


0, 


125 ± 0, 


,014 




4 


0.10 


1 


3.1 


-0.041 




-0.538 




0.105 




0.00 


16 




5 


10.1 ±4.0 


0.121 ±0.048 


-0, 


.441 ± 0.074 


0, 


107 ± 0, 


,037 




8 




4 


7.6 ± 2.8 


-0.013 ±0.028 


-0, 


.527 ±0.057 


0, 


108 ±0, 


,008 


1.00 


4 
16 


0.05 




5 


12.0 ±3.4 


0.014 ±0.016 


-0, 


.009 ± 0.008 


0, 


,109 ±0, 


,027 






0.10 


5 


12.0 ±3.6 


0.007 ±0.035 


-0, 


.009 ± 0.016 


0, 


,101 ±0, 


,026 






0.15 


5 


10.3 ± 1.7 


0.002 ± 0.036 


-0, 


.001 ± 0.020 


0, 


107 ±0, 


,035 




8 


0.10 


1 


4.8 


0.002 




-0.029 




0.097 






4 


0.10 



















0.50 


16 


0.10 


3 


9.8 ± 5.7 


0.009 ± 0.037 


-0, 


,010 ± 0.022 


0, 


130 ±0, 


,009 




8 


0.10 


2 


5.1 ±0.9 


-0.019 ±0.018 


-0, 


,336 ± 0.474 


0, 


124 ±0, 


,001 




4 


0.10 


2 


3.4 ±0.1 


0.003 ± 0.054 


-0, 


,650 ±0.030 


0, 


108 ±0, 


,013 


0.00 


16 




4 


8.8 ±4.0 


0.044 ± 0.088 


-0, 


,507 ±0.088 


0, 


109 ±0, 


,020 




8 




3 


4.1 ±0.9 


0.041 ± 0.039 


-0, 


,588 ± 0.000 


0, 


,114 ±0, 


,004 




4 




2 


5.6 ± 1.9 


0.006 ±0.018 


-0, 


,716 ± 0.289 


0, 


,102 ±0, 


,035 



Note. — (1) Redshift of input group (2) Early-typo fraction of input group (3) Number of input 
galaxies brighter than 7» + 2 (4) Input RS slope (5) Number of mock groups recovered (6) Mean 
number of recovered galaxies in RS (7) Mean recovered RS slope (8) Mean offset of RS effective 
color from input effective color (9) Mean scatter about fit RS 



mg are discussed in iMOfiL The group-finding algorithm 
detects eleven groups in these eight fields, both at zicns 
and elsewhere along the line of sight. Overall, the group 
finding algorithm works exceptionally well, with ten of 
the eleven candidate groups confirmed spectroscopically 
(Table [SJ. The lone exception is the RS at zrs = 0.20 
in the f oreground of the PG1115 lens. While Figure 3 
in IMOItI shows a potential peak at this redshift, there 
are too few velocities to determine if the peak is real. 
We also do not recover two spectroscopic line-of-sight 
groups in our photometry. The spectroscopic group in 
B1422 at z = 0.28 is blended with the lens group RS in 
color space, while the HE2149 group at z = 0.27 simply 
is not detected. 

The photometric redshifts tend to be slightly lower 
than the spectroscopic redshifts (Az = —0.02 ± 0.04), 
because the group finding algorithm tends to return col- 
ors that are slightly bluer than the true RS, as discussed 
in tl3.2l This offset is small, which shows that our photo- 
metric group finding algorithm obtains reliable redshifts 
for detected groups of galaxies. 



Figure ^ shows the spectroscopic group members and 
non- members in comparison with the fit RS for B1422. 
The fit RS agrees with the red edge of the envelope de- 
fined by the group members. Figure [3 shows the spec- 
troscopic group members and fit RSs for both the lens 
group and a background group in PG1115. The back- 
ground group in PG1115 has a significantly higher spec- 
troscopic redshift (zgp = 0.486) than the photometric RS 
redshift (zrs = 0.41), sugg esting that this group may 
have a low /o; indeed. iMOfil find that six of the ten con- 
firmed group members have emission lines. Both of these 
figures illustrate the success of our RS finding technique. 

We define a red-sequence richness parameter TVrs to be 
the number of selected RS galaxies brighter than -I- 2.5 
over and above the normalized background. With two- 
band photometry alone, this is the best richness measure 
we can make. We compare TVrs with the measured ve- 
locity dispersions from iMOfl . Contrary to expectations, 
we find no significant correlation. The explanation may 
be that variable early-type fractions and statistical back- 
ground noise add significant scatter (see also the large 
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TABLE 5 

Parameters of significant candidate red sequences 



RS ID^ 


2R,S 


^spcc 


/. 






(R-I)* 


RS Slope 


Nrs 


RS Scatter 


B0712 : GROUP 1 


0.63 




19.903 






1.361 ± 0.029 


0.115 ±0.029 


9.2 


0.088 


MG0751: GROUP 1 


0.34 


0.3502 


18.229 






0.900 ±0.051 


0.093 ± 0.032 


5.8 


0.044 


MG0751: GROUP 2 


0.48 


0.5605 


19.149 






1.087 ±0.037 


0.121 ± 0.023 


13.1 


0.047 


FBQS0951: GROUP 1 


0.16 




16.461 


1.434 ±0 


060 




0.052 ± 0.042 


5.6 


0.038 


FBQS0951: GROUP 2 


0.27 


0.24'' 


17.675 


1.770 ±0 


041 




0.062 ±0.032 


5.8 


0.047 


FBQS0951: GROUP 3 


0.43 




18.844 


2.349 ± 


037 




0.034 ± 0.024 


6.4 


0.051 


BRI0952: GROUP 1 


0.37 


0.4220 


18.441 






0.938 ± 0.057 


0.059 ± 0.048 


8.4 


0.020 


PG1115: GROUP 3 


0.26 




17.587 


1.732 ±0 


162 




0.063 ±0.081 


4.7 


0.038 


PG1115: GROUP 1 


0.30 


0.3101 


17.924 


1.895 ±0 


045 




0.078 ±0.017 


3.4 


0.031 


PG1115: GROUP 2 


0.41 


0.4859 


18.713 


2.297 ±0 


077 




0.116 ±0.042 


3.4 


0.080 


RXJ1131: GROUP 1 


0.19 


ss 0.104^= 


16.864 


1.508 ±0 


024 




0.059 ±0.020 


25.6 


0.045 


RXJ1131: GROUP 2 


0.29 


0.295*^ 


17.844 


1.852 ±0 


077 




0.071 ±0.046 


6.2 


0.056 


B1422: GROUP 1 


0.30 


0.3387 


17.924 


1.895 ±0 


049 




0.104 ±0.031 


6.0 


0.054 


B1600: GROUP 1 


0.50 




19.262 






1.123 ± 0.064 


0.086 ±0.041 


12.1 


0.035 


MG1654: GROUP 1 


0.22 


0.2527 


17.204 


1.596 ±0 


046 




0.079 ± 0.033 


4.8 


0.021 


B2114: GROUP 1 


0.26 


0.3144 


17.587 






0.792 ± 0.058 


0.055 ±0.030 


7.0 


0.062 


HE2149: GROUP 3 


0.40 


0.4465 


18.647 






0.968 ± 0.060 


0.047 ±0.036 


5.6 


0.052 


HE2149: GROUP 4 


0.59 


0.603 


19.720 






1.288 ± 0.068 


0.086 ± 0.042 


5.3 


0.041 



Note. — The lAU-approved naming convention for these groups is Lens Name: MWKZ GROUP N, where N is the 
group ID; e.g., CLASS BQ712+472: MWKZ GROUP 1 

Matched to'MOf? where appUcable 
''No group data in M06; lens galaxy redshift from Table IT] 

'^Redshift of BRSG from Las Gampanas Redshift Survey IShectman et alj|1996fl 
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Fig. 4. — Comparison of RS-finding algorithm with spectroscopic 
data for B1422. Top panel: Extinction-corrected color-magnitude 
diagram for all galaxies (small crosses) within a projected rvir ra- 
dius surrounding the centroid of the 131422 group. Middle panel: 
Same as top panel, but with best- fitting RS (solid line) and selected 
candidate RS-member galaxies (open diamonds) shown. Bottom 
panel: Same as middle panel, but with spectroscopically confirmed 
group members from IM06i (filled squares, including galaxies from 
beyond Tvir) and confirmed non-members (large crosses) shown. 
This figure shows that the fit RS agrees with the red edge of the 
envelope defined by the spectroscopically-confirmed group mem- 
bers. 



scatter in TV^gee versus a for 2MASS ga laxy clusters with 
cr < 500 km s-i; IKochanek et al.l2003|) . The fuU sample 
of groups from our complete survey will overcome this is- 
sue by spanning a much larger dynamic range, making 
the comparison between N^is and a more meaningful. 



Fig. 5. — Same as FigureBl but for PG1115. Open triangles (cen- 
ter and bottom panels) indicate the candidate RS member galaxies 
for the background group; filled triangles (bottom panel) indicate 
spectroscopic members of the background group. The dashed line 
indicates the fit RS for the background group. This figure illus- 
trates the ability of the red-sequence finding algorithm to detect 
multiple groups in a given field. 

4. RESULTS AND DISCUSSION 
4.1. The Detected Red Sequence Sample 

In total, eighteen significant candidate RSs are de- 
tected in twelve lens fields in the present sample. The 
fit RS parameters are presented in Table [S] and the pho- 
tometry of the selected galaxies is presented in Table IHl 
(available electronically). 

Although some final centroids lie outside rvir of the 
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TABLE 6 

Photometry of candidate red sequence galaxies 



Lens 


Group 


Galaxy 


RA 


Dec 


/ 


SI 


R-I 


S{R - I) 




ID 


ID 














B0712 


1 


8471 


07 16 04.35 


47 08 59.2 


18.61 


0.03 


1.60 


0.03 


B0712 


1 


9253 


07 15 56.86 


47 09 23.7 


20.80 


0.03 


1.45 


0.07 


B0712 


1 


9109 


07 15 58.92 


47 08 52.3 


20.82 


0.03 


1.45 


0.01 


B0712 


1 


9412 


07 15 55.16 


47 09 53.5 


21.21 


0.05 


1.27 


0.11 


B0712 


1 


8204 


07 16 09.09 


47 10 04.5 


21.28 


0.05 


1.10 


0.20 



Note. — (1) Full table is available electronically. (2) The lAU-approved naming 
convention for individual galaxies is: Lens Name: MWKZ GAL NNNN, where NNNN is 
the galaxy ID above; e.g., B0712+472: MWKZ GAL 8471. (3) Units of right ascension 
are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and 
arcseconds {J2000.0). 



lens, we initially detected each RS peak within r^ir of 
the lens galaxy. It is certain that the wide-field images 
contain additional structures outside of this initial search 
radius; no attempt has been made so far to locate such 
structures. 

The average RS slope, A(y-I) /I or A{R-I) /I, for aU 
18 candidate groups is 0.076±0.025. This scatter is com- 
parable to that measured in the mock groups (Tabled, 
so there is no evidence that the RS slope is anything but 
universal at these moderate redshifts. We note that these 
slopes are not corrected for redshift/evolutionary effects. 
Such a correction would require a detailed understanding 
of the origin of the RS and its redshift evolution. 

4.1.1. The Fraction of Lenses in Complex Environments 

Of the twelve lens systems in this paper, eight have 
significant red sequences within Az = 0.06 of the lens 
redshift. Of these systems, six were known or suspected 
previously to have a group at the lens redshift due to 
the need for shear in the lens models or to an overden- 
sity of galaxies on the sky surrounding the lens. The 
two new systems are FBQS0951 (zrs = 0.27) and B2114 
(zrs = 0.26); the B2114 group has been confirmed spec- 
troscopically |M06). 

Figure El shows the CMDs for each of the twelve lens 
fields. In the cases where an RS is detected at the lens 
redshift, the fit RS and the selected RS galaxies are 
shown. Where no RS is detected at the lens redshift, 
a nominal RS for the lens redshift is shown with a slope 
of Pv = 0.09. In many cases, the detected RS is not 
apparent to the unaided eye, but our spectroscopy has 
established that these are physical systems, showing the 
power of our algorithm. 

Figure \7\ presents sky maps of galaxies in each of the 
eight detected RSs near zicns- In many cases, the galaxies 
appear to be clustered, as would be expected in a physical 
group. In a few cases, such as BR10952, t here a ppears 
to be little, if any clustering of galaxies, vet IM06I do find 
a physical group around BR10952 at the lens redshift. 

The sky maps also show that the lens galaxy (indi- 
cated by a star at the origin in each plot) is not nec- 
essarily at the center of the group. Table H presents 
the raw and luminosity-weighted projected centroids for 
each candidate RS, plus (where available) the offset be- 
tween the raw RS centroid and the centroid of spec- 



troscopically confirmed group members from lMQ6t For 
four of the six groups, the centroids agree to within 2a, 
while for the other two (MG0751: MWKZ Group 1 and 
BR10952: MWKZ Group 1), the photometric and spec- 
troscopic centroids disagree. In these cases, the disagree- 
ment could be due to contamination of the photomet- 
ric sample by non-member galaxies, to incompleteness 
in the spectroscopic sample, or to some combination of 
the two. We are analyzing additional spectroscopic data 
for these groups, which should improve the spectroscopic 
completeness and hence the centroid accuracy. 

Among the eight lenses with an RS at the lens redshift, 
the lens galaxy position is within Icr of the projected RS 
centroid for two systems, and within Scr for two more. 
Stated another way, in four cases the lens galaxy is not 
consistent with lying at the centroid (even accounting for 
the centroid uncertainties). This strongly suggests that 
lenses do not necessarily occupy the center of the local 
mass distribution, which is an important realization for 
lens modeling (see ii4.4l andlM06). 

4.1.2. The Fraction of Lenses with Line-of-Sight 
Structures 

Seven of the twelve lens systems have RSs that are 
fore or aft of the lens along the line of sight; prior 
to this study, only one of these groups (B0712) was 
known to have a linc-of-sight group, at z 0.29 
(^assnacht & Lubin 2002). This group is not detected 
here, although we do detect a background group at 
z « 0.67. Two of the hne-of-sight systems (MG0751: 
MWKZ Group 2 and RXJ1131: MWKZ Group 1) have 
very rich RSs indicative of a rich group or even a cluster 
of galaxies, while the remaining six arc more suggestive 
of poor groups. 

Figure |H1 shows the CMDs for all candidate line-of- 
sight groups. We observe a rich, tight RS in RXJ1131, 
indicative of a foreground cluster of galaxies; indeed, 
two of these galaxies have redshift measurements of 
z ~ 0.10 from the Las Campanas Redshift Survey 
()Shectman et al.l I199(t() . and there is extended X-ray 
emission consistent with their position (C. Kochanek, 
private communication). Some of the candidate RSs ap- 
pear dubious, such as the background groups in MG0751 
and HE2149. Once again, however, our spectroscopic 
follow-up confirms the existence of these groups. 

Based on the simulations presented in ii3.2l we expect 
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Fig. 6. — Color-magnitude diagrams and red sequences at the lens redshift. Each panel presents an individual lens. Open diamonds show 
all galaxies within r^ir {= 500h kpc) of the group centroid (or of the lens galaxy for those lenses with no detected red sequence). Selected 
candidate red sequence galaxies are shown as filled diamonds, along with the best-fitting red sequence for galaxies with 7 < /, -j- 2.5 (solid 
line). The lens galaxies are denoted by large, open stars. Asterisks indicate the location of an L» galaxy along the fit red sequence; open 
circles show the location of an L» galaxy at ^lons- For the lenses with no detected red sequence (B0712, B1600, HE2149, and PMN2004), 
a dotted line indicates a nominal red sequence with a slope of 0.09 at zicns- Eight of the twelve lenses have a red sequence at the lens 
redshift. 



about one false positive (either line-of-sight or at zjcns) 
among the five lens fields observed in the l/-band and 
two to three false positives from the seven lens fields with 
i?-band data. 

4.2. Lens Galaxy Environments and Shear 

It is interesting to consider whether observed lens 
galaxy environments can explain the external tidal shears 
required in lens models. Shear is created by an asym- 
metric distribution of mass around the lens galaxy; a 
simple way to get that is to have the lens galaxy off- 



set from the centroid of the group. Historically, several 
attempts to compare the shear required by lens models 
with the distribution of galaxies within ~ 140 kpc of 
the lens have yielded mixed results. In two quadruple- 
image lenses (PG1115 and B1422), the model shear ap- 
pears to be co nsistent with the distribution of galaxies 
near t he lens llHogg k Blandfordlll994t iSchechter etai\ 
19971 IKeeton fc Kochanekl 119971 TKundic et all Il997al: 
Kundic et alJll997bj) . In contrast, in the Einstein ring 
MG0751, the shear angle i s not consistent w ith the im- 
mediate lens environment l|Lehar et alJll99'^ Also, in a 
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sample of 10 double-image lenses. iLehar et all 1)20001) do 
not find any strong correlation between model shears and 
close environments. 

We can improve upon the previous studies in three 
ways. First, we have obtained a larger sample of lens 
environments. Second, we have measured lens environ- 
ments over a larger area on the sky. Third, we have 



identified enough group member galaxies to make a re- 
liable measurement of the group centroid position. One 
limitation of our present analysis is that, without mea- 
sured velocity dispersions, we cannot estimate the am- 
plitude of the shear contributed by each lens environ- 
ment. (See M06, for an analysis of shear amplitudes for 
a smaller sample of lenses.) Nevertheless, we can still 
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Fig. 7. — Sky maps of red sequence members of candidate lens groups. Point sizes are scaled by galaxy magnitude relative to /«. The 
open star at each map's origin marks the lens position, and the error bars denote the positions of the unweighted red sequence centroids 
and their Itr errors. The circle shows one projected group virial radius (= 500/i~^ kpc) at the red sequence redshift centered on the group 
centroid. These maps show that, in many cases, the lens galaxy is not consistent with the group centroid and is not the brightest galaxy 
in the red sequence. 



compare the shear position angle with the position angle 
of the group centroid (relative to the lens galaxy) and de- 
termine whether the observed environments are at least 
consistent with the model shears. 

TablelHlgives the position angles for the observed group 
centroid and the lens model shear, where available. In all 
six cases where both are available, the centroid and shear 
position angles are remarkably consistent (albeit with rel- 
atively large centroid errorbars in a few cases). Perhaps 
the most interesting system is MG0751: while the imme- 
diate lens environ ment cannot explain the model shear 
ijLehar et al.lll99'^ . we now see that the larger environ- 
ment can. This example illustrates why it is important 
to observe lens environments over an area large enough 
to include the full virial extent of a group at the lens 
redshift. It will be interesting to reconsider this ques- 



tion in more detail, after measuring velocity dispersions 
so that we may estimate shear amplitudes as well as po- 
sition angles, and after making new lens models of all 
the systems. Still, it is very encouraging to see new ev- 
idence that the observed lens environments can explain 
the shears required by lens models. 

4.3. Lens Galaxy Magnitude Distribution 

The sample of lenses allows us to examine the dis- 
tribution of lens galaxy magnitudes. While bright 
(i.e., massive) galaxies have large lensing cross-sections, 
there are many more faint galaxie s in groups (e.g., 
iKeeton. Christlein. fc Zabludofll2000|) . Taken together, 
these effects lead one to ex pect t he typical lens to be 
~ ijKochanek et al.1 12000'. 20011. though it is stiU un- 
clear how faint (and dwarf-likc) lenses can be. In our 
sample, ten of the eleven lens galaxies with redshift mea- 
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surements are fainter than (see Table E))- In order to 
determine whether our sample includes some unknown 
bias toward faint lenses (although one would expect that 
any bias would be toward bright, massive galaxies), we 
calculate the expected distribution of lens galaxy mag- 
nitudes based on theoretical expectations and a minimal 
number of reasonable assumptions. We then compare 
this theoretical distribution with the observed distribu- 
tion in our sampl e and a larger sam ple of gravitational 
lens galaxies from lRusin et al.l (|2003). 

Lensing's sensitivity to mass means that the distribu- 
tion of lens galaxy magnitudes is not the same as the 
galaxy luminosity function. Nevertheless, it is straight- 
forward to pr edict the magnitude distr ibution for lens 
galaxies fe.g.. iFukueita fc Turner! [19 91*). For a singu- 
lar isothermal sphere lens model, the lensing probability 
F oc gf , where aj, is the velocity d ispersion of the lensing 
galaxy l)Fukugita &: Turner! 19911) . A ssumin g a Fa ber- 
Jackson relation L oc a'*^-' l|Faber fc Jackson. ll976j) , we 
derive a lensing probability oc L^/tf' . 



We beg in with a Schechter luminosity function 



We beg m wit 
llSchechteijll97l . 

ML) dL 



(10) 



Factoring in the lensing probability, the probability for a 
lens galaxy to have a luminosity L is 

gi[L)AL = ^(^] e-^/^-di. (11) 



We integrate this equation and normalize to a limit- 
ing luminosity Lunj (thereby accounting for bias due 
to magnitude-limited surveys and non-converging func- 
tions) to get the cumulative probability for lens galaxy 
luminosities greater than i/L*, 

^ ' l-r(H-a + 4/7Fj,Liim/i*) ^ ' 

where F is the normalized incomplete gamma function. 
This equation is readily modified to / magnitudes using 
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TABLE 7 

Centroids of candidate red sequences 



RS ID Raw Centroid L- weighted Centroid Offset from M06 









o 




I 








a 










\ ) 


B0712: GROUP 1 




07 


16 


03 


47 


09 


25 


07 


16 


03 


47 


09 


14 




MG0751: GROUP 1 




07 


51 


48 


27 


16 


42 


07 


51 


47 


27 


16 


45 


101 ± 20 


MG0751: GROUP 2 




07 


51 


44 


27 


17 


10 


07 


51 


44 


27 


17 


09 




FBQS0951: GROUP 


1 


09 


51 


22 


26 


35 


01 


09 


51 


23 


26 


35 


10 




FBQS0951: GROUP 


2 


09 


51 


20 


26 


36 


19 


09 


51 


19 


26 


36 


19 




FBQS0951: GROUP 


3 


09 


51 


22 


26 


36 


01 


09 


51 


23 


26 


35 


51 




BRI0952: GROUP 1 




09 


55 


03 


-01 


30 


02 


09 


55 


03 


-01 


30 


01 


107 ± 41 


PG1115: GROUP 3 




11 


18 


16 


07 


43 


27 


11 


18 


16 


07 


43 


12 




PG1115: GROUP 1 




11 


18 


18 


07 


45 


27 


11 


18 


17 


07 


45 


47 


20 ±22 


PG1115: GROUP 2 




11 


18 


18 


07 


46 


03 


11 


18 


19 


07 


46 


24 




RXJ1131: GROUP 1 




11 


32 


00 


-12 


31 


48 


11 


32 


03 


-12 


32 


03 




RXJ1131: GROUP 2 




11 


31 


50 


-12 


31 


45 


11 


31 


49 


-12 


31 


45 




B1422: GROUP 1 




14 


24 


38 


22 


55 


54 


14 


24 


37 


22 


56 


15 


40 ±24 


B1600: GROUP 1 




16 


01 


45 


43 


16 


48 


16 


01 


45 


43 


16 


44 




MG1654: GROUP 1 




16 


54 


37 


13 


46 


38 


16 


54 


40 


13 


46 


39 


54 ±39 


B2114: GROUP 1 




21 


16 


47 


02 


24 


54 


21 


16 


47 


02 


24 


41 


68 ± 40 


HE2149: GROUP 3 




21 


52 


08 


-27 


32 


17 


21 


52 


06 


-27 


32 


09 




HE2149: GROUP 4 




21 


52 


06 


-27 


31 


31 


21 


52 


05 


-27 


31 


36 





Note. — Units of right ascension are hours, minutes, and seconds, and units of declination 
are degrees, arcminutes, and arcseconds (J2000.0). 



TABLE 8 

Comparison of group centroids to lens properties 



Lens 


Raw Centroid 


L-weighted Centroid 


Offset 


Position Angle 


Shear 
PA 


Reference 




Aa 


AS 


Aa 


AS 


(kpc) 


Raw Centroid 


L-weighted 


MG0751 


90± 18 


10 ±15 


76 ± 26 


13± 14 


440 ± 90 


84± 10 


80± 14 


64 


(1) 


FBQS0951 


-37± 22 


65 ±20 


-53 ± 32 


65± 16 


310 ±90 


-30 ±17 


-39 ±19 




BRI0952 


47± 21 


4± 16 


48 ± 18 


6± 18 


240 ± 110 


85 ±27 


83 ±24 


66 


(2) 


PG1115 


15± 16 


-32 ± 27 


-4 ±22 


-12 ±20 


160 ± 110 


155 ± 47 


-162 ± 78 


-113 


(3) 


RXJ1131 


-21 ± 23 


14 ±20 


-35 ± 21 


13± 15 


110 ± 100 


-56 ± 58 


-70 ±35 


-80 


(4) 


B1422 


2 ±22 


-7± 17 


-17± 18 


15± 17 


30 ±80 


164 ± 86 


-49 ± 54 


-53 


(4) 


MG1654 


-76 ±41 


16 ± 16 


-23 ± 37 


17± 18 


280 ± 140 


-78 ±28 


-53 ± 63 


-80 


(5) 


B2114 


-58± 17 


-53 ±21 


-52 ± 18 


-66 ± 22 


320 ± 80 


-132 ± 15 


-142 ± 14 





References. — (1) ILehar et alHT99^ (2) ILehar et al1l2000l (3) llmpev et aLllTggl (4) IKeeton. Gaudi. He Pett'er3l200l (5) 
iKochanck 1995 

Note. — Units of Aa and A<5 are arcseconds (J2000.0). PA measured East of North. Except for PG1115, shear PAs are 
defined modulo 180 degrees. 



L/L^ ^ IQ-O-m-i,)^ Note that Eq.Ulis fully degenerate 
in a and 7fj. 

We compare two samples to this theoretical prediction. 
First, we use the nine lenses from the present sample that 
have early-type lens galaxies (we exclude the two spi- 
ral lenses B1600 and PMN2004 and the two-plane lens 
B2114). Second, we note that as this simple theoretical 
prediction does not explicitly involve environment, we 
can also compare it wit h the larger s ample of early-type 
lens galaxies defined bv lRusin et alJ 1^2003.) . who do not 
determine the environments of their lenses. The Rusin 
sample spans a larger redshift range than ours, so we use 
only the low-redshift systems (zions < 0.6) with spectro- 
scopic redshifts, a total of 12 lenses. Five of the nine 
lenses in our sample overlap with the Rusin sample. 

We take obser v ed F 814W or F791W magnitudes 
from iRusin et all pOO^ or from the CASTLES web 
database,^ convert to Cousins /, correct for Galactic ex- 

^ In most cases the lens galaxy magnitude cannot be determined 
from our data, because the lens galaxy and quasar images are un- 



TABLE 9 
Lens Galaxy Magnitudes 



Lens 


-^Icns 


-fo.lcns — I* 


B0712 


19.58 ± 0.07 


0.64 


MG0751 


21.28 ± 0.64 


3.05 


FBQS0951 


19.68 ± 0.23 


2.01 


BRI0952 


21.23 ± 0.04 


2.79 


PG1115 


18.94 ±0.02 


1.02 


RXJ1131 


16.6 ±0.3 


-1.3 


B1422 


19.68 ± 0.25 


1.76 


B1600 


20.80 ±0.41 


1.54 


MG1654 


17.92 ± 0.02 


0.72 


PMN2004 


20.51 ± 0.02 




B2114 


18.65 ± 0.23 


1.06 


HE2149 


19.60 ±0.03 


0.87 



Note. — Except for PMN2004, all 
lens galaxy magnitudes are from the lit- 
erature (see text). 
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Fig. 8. — Color-magnitude diagrams for candidate line-of-sight red sequences. Symbols are as in Fig. |^ Ten line-of-sight lenses are 
found in seven (of twelve total) lens fields. Although to the naked eye many of these red sequences are not obvious, our spectroscopic 
observations of HE2149, MG0751, and PG1115 confirm the line-of-sight groups in HE2149 and MG0751 and the background group in 
PG1115 (Momcheva et al. 2006). The low-z line-of-sight candidate in PG1115 was not detected spectroscopically. The apparently steeper 
red-sequence slopes for the lower four panels are due to the smaller vertical color range in R — I. 



tinction, and calculate I — I^. For we again use our 
photometric models of early-type galaxies evaluated at 
the RS photometric redshift (if an RS exists) or at the 
lens spectroscopic redshift. For our lens sample, we use 
the photometric RS redshifts to insure consistency with 
our ongoing, larger lens sample, in which many lens spec- 
troscopic redshifts will not be known. The use of RS 
photometric redshifts instead of spectroscopic redshifts 

resolved. 



does not affect any of our conclusions here. 

The cumulative distributions of observed lens galaxy 
magnitudes are shown in Figure El along with predicted 
distributions for different values of a and 7fj. Be- 
cause our sample extends fainter than the Rusin sample 
{Iiim — « 3 versus 1.5, respectively), the two observed 
samples are plotted separately. There is, however, no 
significant difference between our sample and the Rusin 
low-redshift subsample: the Kolmogorov-Smirnov (KS) 
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-1 



Fig. 9. — Cumulative fraction of lens gal axy magnitudes com- 
pared with predicted distributions from Eon. I12l Solid curves rep- 
re sent 7f,t = 4 ; dash ed curves arc for 7pj = 3.3, as calculated 
in IRusin et al.l 120031) . From left to right, curves are for faint- 
end luminosity function slopes of a = —0.6 to — 1.8 in intervals 
of 0.3, with the middle value of —1.2 shown in bold. Top panel: 
Lens galaxies in this paper's sample (solid points), with theoreti- 
cal curves nor malized to I — /, =3. Bottom panel: zions ^ 0.6 
lens sample of lRusin et al.l |2003) (open circles). The lens galaxy 
magnitude distribution in both samples is well-described by the 
theoretical models and reasonable values of 7pj and a. 



probability is P = 0.42. Both samples indicate that the 
median low-redshift, early- type lens galaxy has a mag- 
nitude / — /* ss 0.75, or a luminosity L « 0.5L*. This 
median value is directly dependent on the value of /*; 
as discussed in Appendix A, our normalization of /, is 
brighter by 1 magnitude than most modern surveys, 
so the median lens luminosity may well be closer to L,. 
Even allowing for this, fully half of the elliptical lens 
galaxies are fainter than L», some significantly so (in 
particular, BRI0952 and MG0751 are 2 to 3 magnitudes 
fainter than having L ~ O.IL,). 

The figure shows that there are values of a-|-4/7Fj lead- 
ing to theoretical predictions that match the data reason- 
ably well. High values of a -I- 4/7fj are clearly inconsis- 
tent with the data: KS tests rule out a + 4/7fj > — 0.2 
at the 2(7 level (P < 0.05). For reasonable Faber- Jackson 
slopes (7FJ ~ 3 to 4), this would suggest that the faint- 
end luminosity function slope is steep, a < —1.2. We 
note that any selection bias against small image separa- 
tion lenses will tend to lessen the slope of the observed 
distribution in Figure 9, as such a bias will be more pro- 
nounced in the less massive (and, by assumption, fainter) 
lenses. This potential bias, the degeneracy between a 
and 7FJ, and the uncertainty in the normalization in /, 
prevent us from drawing detailed conclusions here. 

There are several important results from this exercise. 
On the observational side, the median lens galaxy lu- 
minosity is near or below i*, and the luminosity distri- 
bution is very broad. In other words, lens galaxies are 
not all massive, super-L* ellipticals; rather, they span 
a range of luminosities (and, presumably, masses), ex- 
tending down at least to ~ O.IL*. On the theoretical 
side, the observed luminosity distribution is consistent 



with simple models of a Schechter luminosity function 
weighted by the lensing probability, for reasonable val- 
ues of a and 7fj. The key conceptual point is that, while 
lensing selects galaxies by mass, dwarf galaxies are suf- 
ficiently common to comprise a moderate fraction of the 
lens galaxy sample. 

4.4. Properties of the Brightest Red Sequence Galaxies 

In the local universe. X-ray luminous groups 
of galaxies with a > 300 km s^^ contain a gi- 
ant elliptical, the Briqhtest Group Galaxy or BGG 
(^ulchacv & Zabludofi 1998), which occupies a unique 
position in the kinemati c and projected s patial 
center of the group ( . Zabludoff fc MulchaevI Il998t 
Ivan den Bosch et all 1200,'tI) ^ In this section we exam- 
ine the BGGs in our group sample in order to deter- 
mine whether any evolution in their properties can be 
observed. In two of three distant, high-cr (~300-500 km 
s~^) groups, IM06. found that the lack of a significant 
offset of the BGG from the group centroid is similar to 
that observed in nearby groups. Given the larger sample 
offered by this photometric study, we identify the BGG 
in each candidate group and examine its relation to the 
group centroid. As BGGs at low redshift are typically el- 
lipticals, we assume here that the brightest galaxy on the 
selected RS (BRSG, for brightest red sequence galaxy) 
is the BGG. 

The properties of the BRSG in each group are given in 
Tablelini Only two lens galaxies, RXJ1131 and MG1654, 
are BRSGs. We also note that two BRSGs, those in 
B2114 Group 1 and HE2149 Group 4, are known not 
to be the BGG as their redshifts are inconsistent with 
group membership. The fraction of lens galaxies that 
are B RSGs is 0.25 (two lenses out of eight lens groups); 
iMOfil find a BGG fraction of 0.33 (two lenses out of six 
lens groups). These values compare favorably with the 
suggestion by Oguri ( 2006) than < 50% of lens galaxies 
occupy the dominant halo, the remainder being "satel- 
lite" galaxies, or sub-halos within the overall group halo. 
This result is also not surprising in light of our conclu- 
sion in i)4.3l that the lens galaxy is typically not a very 
luminous galaxy. 

Figure EH compares the offset of each BRSG from 
the raw projected spatial RS centroid to the measured 
group velocity dispersion from lMO^ With one exception 
(MG0751), groups with a > 300 km s"^ have a BRSG 
near the group centroid, while none of the smaller-a 
groups have a central BRSG. A Spearman rank corre- 
lation of a and the centroid-BRSG offset also suggests a 
trend (the probability of no correlation is 0.14). Our data 
therefore intimate that there are dynamically evolved 
(i.e., virialized) groups with high a and a central, domi- 
nant early-type member by z ~ 0.5, groups whose prop- 
erties are consistent with evolved (virialized) groups in 
the nearby Universe. 

4.5. The Ratio of Quad to Double Lenses 

One unsolved puzzle in lensing studies is why there 
are so many quadruple-image lenses compared with 
double-image lense s (Kine ct al. 1996; Kochanek 19963 
Kecton & Kochanek ' 1997; Rusin & Tcemark 20o3 
Cohn fc Kochanek ,2004: Keeton fc Zabludoff 2004i) . 
Highly-flattened lens galaxies are one way to obtain 
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TABLE 10 

Properties of brightest red sequence galaxies 



RS ID 


IDbrsg 


^o.brsg 


Offset 


Tom lens 


Offset from ccntroid 


IDbgg" 


Notes 








(") 


(kpc) 


("1 








T3n'710. /^"DOTT"D 1 


8471 


18.387 


■ ■ ■ 


■ 


28 ± 10 


190 ± 70 


■ 




IVitjU/oi; CjHUU-r 1 


8626 


17.824 


6 


30 


96 ± 18 


460 ± 90 


8626 


a.K.a. Lri ( ionrv & tvocnaneK lyyy) 




8387 


18.171 






34± 12 


200 ± 70 


8476 


rJrignter tnan^^Ub -dLtCj 


T?T3/~\Gnnc;i . r^T^OTTr) i 


( ooy 


1 7 "71 Q 


■ ■ ■ 


■ ■ ■ 


128 ± 49 


350 ± 130 






rrito2oOyol: CjKUU-T 2 


9345 


16.941 


110 


450 


65 ± 22 


270 ± 90 






T?TI?r~\Cnnc; 1 . r^T50TTT3 Q 

rrS^bOyol; CjKUU-t o 


8679 


17.665 






18 ± 20 


100 ± 110 






BRI0952: GROUP 1 


9622 


19.422 


71 


360 


84 ± 18 


430 ± 90 


8406 


Brighter than M06 BGG 


PG1115: GROUP 3 


11679 


18.804 






106 ± 21 


430 ± 90 






PG1115: GROUP 1 
PG1115: GROUP 2 


12600 
13804 


17.824 
17.732 


24 


110 


42 ± 19 
101 ± 17 


190 ± 90 
550 ± 90 


12600 
13764 


a.k.a. Gl flmpey et ah 1998) 
Brighter than M06 BGG 


RXJ1131: GROUP 1 


7236 


15.135 






44 ± 22 


140 ± 70 




z = 0.104 from LCRS 


RXJ1131: GROUP 2 


6323 


16.6 








25 ± 22 


110 ± 100 




Lens galaxy 


B1422: GROUP 1 


13152 


17.523 


8 


40 


2 ± 22 


10 ± 100 


13152 




B1600: GROUP 1 


7950 


19.714 






88 ± 15 


540 ± 90 






MG1654: GROUP 1 


17334 


17.92 








78 ±40 


280 ± 140 


17334 


Lens galaxy 


B2114: GROUP 1 


7943 


17.841 


185 


740 


108 ± 20 


430 ± 80 


8938 


Not group member ^0^ 


HE2149: GROUP 3 


7568 


17.663 






80 ± 15 


430 ± 80 


8364 


Brighter than M06 BGG 


HE2149: GROUP 4 


7534 


19.892 






54± 16 


360 ± 110 


6342 


Not group member (MQ^ 



^ID of BGG inlMOi 



200 400 600 

"Group (km/s) 

Fig. 10. — Offset of the brightest red-sequence galaxy (BRSG) 
from the red sequence centroid as a function of spectroscopic group 
velocity dispersion (fromlMOg). Crosses represe nt B RSGs known 
not to be the brightest group galaxy (see Table . The dotted 
box indicates the maximum spatial offset of BGGs in groups with 
O"e-r oun > 300kms~^ in the local Universe (Zabludoff & Mulchacv 
1199 81. As in the local Universe, higher velocity dispersion groups 
have a BGG near the group centroid, while lower velocity disper- 
sion groups groups do not. 



groups around the lens galaxies. 

Photometrically, we find some weak evidence that quad 
lenses may have more complex environments than dou- 
ble lenses. Three of the four quad lenses and all three 
"other" lenses have associated RSs, while only two of the 
five double-image lenses have RSs. However, this differ- 
ence is not statistically significant; larger samples are 
needed to make any definitive statement. The lack of 
strong evidence for a difference between quad and dou- 
ble lens environments contrasts with our findings in iM06l 
The difference may simply reflect small number statis- 
tics. A more likely explanation is that the results in M06 
were based on detailed calculations of the convergence 
and shear due to each lens environment — calculations 
we cannot make in this paper because we do not know 
the true memberships and velocity dispersions for all of 
our groups. Our complete lens sample will enable us to 
reexamine any connections between image number and 
environment that may be hidden in the noise of our cur- 
rent sample. 

We find no correlation between the number of lensed 
images and the presence/absence of line-of-sight struc- 
tures. Three of the four quads, three of the five doubles, 
and two of the three "others" have at least one line-of- 
sight structure identified photometrically. Again we need 
more data to draw a strong conclusion about any possible 
relation. 



the seemingly high numbers of quad lenses; another 
possibility is the presence of a large tidal shear from 
mass near the lens galaxy or along the line-of-sight. In 
this latter scenario, quad lenses might be expected to 
be located preferentially in complex environments or 
along cluttered lines of sight. Indeed, IM06I note that 
in a sample of eight lenses, both of the quad lenses 
and the one quad/ring lens are in relatively massive 
groups, while another ring and two of the doubles are 
in poorer groups, and the other two doubles have no 



5. CONCLUSIONS 

We have presented the first results from an on- 
going wide-field imaging survey of ~ 80 strong gravi- 
tation al lens sy stems. We use a modified version of the 
■Gladders fc Yea (,,2000.) group-finding algorithm to locate 
potential red sequences, both at the lens redshift and 
along the line-of-sight, in the fields of twelve lenses. This 
modified algorithm is capable of detecting a > 200 — 300 
km s~^ poor groups out to z « 0.4, and higher-u groups 
out to higher redshifts with a low rate of false positives 
and reasonable completeness. 
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We draw the following conclusions: 

• Most gravitational lenses lie in complex environ- 
ments. For this sample, 67% of lens galaxies (eight 
of twelve) lie in galaxy groups. As there are no 
strong biases in this sample, is is unlikely this frac- 
tion would drop below ~ 50% once the environ- 
ments of all lenses are surveyed. 

• Most gravitational lenses have interloping struc- 
tures. In particular, we detect ten groups projected 
within ^ 1.'5 of seven lenses; i.e., « 60% of the 
lenses in this sample have at least one interloping 
group that could impact the lens model (see also 
Momcheva et al. 2006). 

• The centroid positions of the lens groups are con- 
sistent with the directions of the external shears 
required by lens models, for the six systems where 
we can make the comparison. This suggests that 
lens environments can explain the model shears — 
which has often been assumed, but not well proven. 
In at least one system (MG0751), the agreement is 
seen only when the full virial extent of the group is 
included. 

• The typical gravitational lens is not a super-L* el- 
liptical. Rather, the observed distribution of lens 
magnitudes is well-describ ed by a convolut ion of 
lensing probabilities with a lSchechted l)1976f ) lumi- 
nosity function, with a median lens galaxy lumi- 
nosity ~ and with known lens galaxies as faint 
as - O.IL*. 

• The typical lens galaxy is not the brightest group 
galaxy. In the current sample, only two of eight 
lens galaxies (25%) in a group are the brightest 
group galaxy. 

• As in the local Universe, higher velocity dispersion 
(cr > 300 km s""'^), intermediate-redshift galaxy 
groups have a brightest group galaxy near the 
group centroid, whereas the brightest group galaxy 
typically lies outside the group center in lower ve- 
locity dispersion groups. This result suggests that 
the higher a groups are more dynamically evolved 
(i.e., viriaHzed) than the lower-cr groups and that 
evolved (virialized) groups exist by z 0.5. 

• In total, we report the photometric discovery or 
recovery of 18 candidate poor groups of galaxies 



in the redshift range z ^ 0.2 to 0.6. The average 
red sequence slope is 0.076 ± 0.025, in agreement 
with that of nearby groups. The measured scat- 
ter about the red sequence is comparable to that 
expected from random photometric errors. When 
analysis of the full lens sample is complete, we ex- 
pect to find ~100 groups at these intermediate red- 
shifts, permitting detailed studies of the evolution 
of group structure and members. 
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APPENDIX 



PREDICTED GALAXY PHOTOMETRIC PROPERTIES 

For this paper, we determine the photometric properties of galaxies as a function of redshift from simple evolutionary 
models alr eady used for analy s is of l ens galaxies (Kccton ct al. 1998). These models begin with spectrophotometric 
models bv iBruzual fc Charloti l)1993|) . assuming solar metallicity and a star formation epoch sta rting at z — 5. We 
assume that early- type galaxies are formed by a 1-Gyr burst of star formation with a ISalpet'e^ (t l955) i nitial mass 
function (IMF) followed by passive evol ution to the present day. We mode l late-t ype galaxies with a lScalol ljl986j) IMF 
and with the star formation rates from iGuiderdoni fc Rocca-Volmerangd lll988D for Sa, Sb, and Sc galaxies. We use 
these models to compute the color-, and evolutionary corrections needed to predict the apparent magnitude and 
colors of an galaxy as a function of redshift. For this work, we define to be the luminosity of a z = galaxy 
with an absolute i3-band magnitude of Mb(L*) = —19.9 -I- 51og/i. The model color-magnitude curves are shown in 
Figure [TTl 
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Fig. 11. — Model galaxy V — I (left) and R — I (right) colors versus I magnitude for an L» galaxy. Curves indicate E/SO (solid), Sa 
(dotted), Sb (dashed), and Sc (dot-dashed) galaxies and range from z = 0.1 — 1.0. Boxes mark redshift intervals of Az = 0.15, with the 
corresponding numerical redshift indicated along the E/SO track. 



We also consider the impact of different normalizations of Mb{L^,) on our results. We have converted each of the 
published normalizatio ns to a consistent cosmology and to Mr , the latter by using the z — colors for elliptical 
galaxies calculated by Fukugita. Shimasaku. & Ichikawa' 11995j) . Recent analysis of the density-dependence of the 
galaxy luminosity function from the 2dF Galaxy Redshift Survey (CoUcss et al. 200]J) is presented in Croto n et alJ 
()2005j) . who find Mb{L^) ranging from —18.3 -f 5 log ft, for early- type galaxies in voids to —1 9 .8-1- 5 1og/i for early- 
type galaxies in clusters, with a value ~ —19.4 -t- blogh for the entire volume. iBernardi et al.l l)2005() find a value of 
Mb{L^) = — 18.9 + 51og/i for early-type galaxies from the Sloan Digital Sky Survey. Other galaxy surveys have values 
of Mb{L^) ranging from w -18.7 -|- 5 log /i to « -19.4 5 log ft (e.g., see Table 1 in lGonzalez et al.ll2000j) . 

Analytically, the change in the color V — I oi a galaxy of magnitude / on the red sequence from a change in 
normalization of of A/» is: 

A(y - /) = l3vAh , (Al) 

where /3y is the slope of the red sequence mV — I. Taking the faintest value of Mb{L:,) — —18.3 + 5 log ft from the 
surveys above and (3v = 0.08 (30J, the maximum change in the color of a red sequence galaxy is A{V — I) = —0.128. 
Over the redshift range of 0.2 < z < 1.0, the average change in photometric redshift as a function of change in color is 
Az/ A{V — I) = 0.133, so the maximum change in redshift from a change in normalization would be Az ~ —0.02. We 
test this calculation by running the group-finding algorithm on several groups with AlBiLsr) = —18.3 + 5 log ft, and 
indeed the detected group redshifts are lower by ~ 0.02. We therefore note that our photometric group redshifts may 
be systematically biased toward higher redshifts by as much as Az = 0.02, smaller than the measured scatter between 
groups photometric and spectroscopic redshifts given in ti3.3l 
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